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Netermination of Carbon 14 in the Terminal Positions of 
Sugars: Preparation of p-Arabinose-5-C* From 


v-Fructose-1,6-C* “” 


Harriet L. Frush and Horace S. Isbell 


Methods are described for the determination of carbon 14 in the terminal positions of 


reducing sugars 


Phe sugar in alkaline solution is oxidized by 


means of molecular oxygen 


to the next lower aldonic acid, and the acid is separated in the form of a convenient deriva 


tive. 


activity located at this carbon. 


Phe activity 


The loss in radioactivity coincident with the removal of carbon | is a measure of the 
of the carbon at the nonreducing end of the 


molecule is determined by periodate oxidation of the aldonic acid, followed by formation and 
radioactivity assay of the dimedone compound of the resulting formaldehyde 


The method was applied to the labeled p-fructose prepared from 
suboridans, and the p-araboniec acid formed from the p-fructose 


by the action of Acetobacte 
was isolated as potassiura D-arabonate 


It was found that all of the 


p-mannitol-1-C' 


f activity of the labeled 


p-fructose was located equally at carbons 1 and 6 


Potassium p-arabonate-5-C"™ was 
26-percent radiochemical yield 
‘ 
it 


\ ield 


p-fructose-1,6-( was ther 


vield of 18 percent of the original p-mannitol-1-¢ 


1. Introduction 


Previous reports from this laboratory |1, 2] ° have 
escribed the preparation of 1-C'-labeled mannose, 
e reduction of this to p-mannitol, and the oxida- 
on of the p-mannitol by A, suboridans to labeled 
fructose. As the two halves of the mannitol 
wlecule are stereomerically alike, oxidation by 
.suboridans can take place either at carbon 2 or at 
arbon 5 If oxidation took place without either 
n isotope effect or breakdown of the p-mannitol and 
synthesis of the fragments, the derived p-fructose 
ould be labeled equally and exclusively at carbons | 
nd 6. The present report presents evidence that 
his fructose is, within experimental error, labeled 
jually and exclusively at carbons 1 and 6, and de- 
ribes convenient methods for determining carbon 
4 at these positions. 


2. Discussion of Experimental Work 


Ordinarily, considerable difficulty is experienced 
i the purification of p-fructose. In the course of 
ie work on the labeled sugar, a convenient proced- 
re was developed for the crystallization of p-fructose 
rom crude sirups. This consists of dehydration of 
he sugar sirup by azeotropic distillation with abso- 
te ethanol, followed by dissolution in about two 
arts of methanol, and addition of a higher alcohol 
ich as butanol or isoamyl alcohol to the point of 
heipient turbidity. After nucleation, the solution 
lowly deposits large crystals of p-fructose. 
on the development of methods for t 
ponsored by the Atomic Energy C 
fructose-1,6-C'™”’ covers all materials labeled 


proportion of the labels As the f 
mannitol-1-C™, positions 1 and 6 are not 


wckets indicate the literature reference 


also prepared from pb-mannitol-1-¢ 
without the intermediate separation of 
converted to p-arabinose-5-C'* in an 


‘ou 52-percent 
ervstailine 
over-all radiochemical 
‘ 


For assay of carbon 1 of p-fructose, the method 
chosen was an alkaline degradation of the sugar in 
the presence of molecular oxygen |8, 4] adapted to a 
semimicro scale and so conducted as to avoid con- 
tact of the sugar and alkali except under oxygen 
pressure. The resulting p-arabonic acid was con- 
veniently separated as crystalline potassium p-ara- 
bonate, and the assay of carbon 1 of the labeled 
p-fructose was obtained fron the difference in activity 
of the sugar and the potassium salt. The method 
is convenient, and can be used for the assay of car- 
bon 1 of any simple aldose or 2-ketose by choice of 
a suitable derivative for separation of the lower 
sugar acid. For the assay of carbon 6 of p-fructose, 
the potassium p-arabonate described above was oxi- 
dized with sodium periodate, and the formaldehyde, 
which was derived only from the terminal carbinol 
carbon, was isolated directly from the oxidation 
mixture by means of the dimedone compound. The 
procedure is somewhat similar to that used by 
Reeves [5] for the quantitative determination of the 
primary carbinol groups of sugars; it does not require 
distillation of formaldehyde prior to formation of 
the dimedone compound. The method is suitable 
for the assay of the terminal carbinol group of any 
aldonic acid or simple sugar. Although aldoses can 
be assayed directly by the periodate oxidation 
method described here, direct assay of ketoses is not 
practicable because they contain two primary car- 
binol groups. Thus assay of ketoses by alkaline 
oxidation to the lower aldonic acid, and periodate 
oxidation of this, is particularly advantageous. 

The C-labeled p-frictose, assayed as described 
above, yielded potassium p-arabonate with a mole- 
cular loss of 50 percent of the activity. The for- 
maldehyde-dimedone derivative obtained after peri- 
odate oxidation accounted for all of the activity of 
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Thus the original 


p-arabonate 
and exclusively in 


the potassium 
p-fructose was labeled equally 
positions | and 6 

The methods for the assay of carbon | of p-fructose 
ean be used without change for the assay of glucose 
and mannose. By alkaline degradation in the pres- 
ence of oxygen, it was shown that all of the activity 


of a sample of p-glucose-1-C" prepared by the 


evanohydrin synthesis from p-arabinose and radio- 
active cyanide is, unquestionably, located at carbon 1. 


The semimicro procedure developed for the prep- 
p-arabonate-1-C™ is useful, 


aration of potassium 
p-arabonic-5-C™ acid 


also, for the preparation of 
and p-arabinose-5-C™. In the preparation of these 
s unnecessary to separate crystalline 


substan es it} 
sub- 


p-fructose after oxidation of p-mannitol by A. 
oridans The fermentation mixture, after depro- 
teinization, deionization, and alkaline oxidation, 
gave potassium p-arabonate in 52.4-percent yield, 
based on the p-mannitol fermented However, be- 
cause one-half of the activity Was lost in removal of 
the carbon 14 recovered as potassium 
was 26.2. percent. The potassium 
was converted to p-arabonic-5-C"™ 


carbon 1, 
p-arabonate 
p-arabonate-5-C"™ 
lactone and this was reduced to 
in 69-percent yield by the procedure developed for 
the preparation of p-arabinose-1-C" [6]. Obvi- 
ously, the 5-labeled p-arabinose provides a starting 
material for the preparation of p-glucose-6-C™ and 


- ae 
p-arabinose-5-C" 


p-mannose-6-C 


3. Experimental Details 
3.1. Purification of p-Fructose 


For purification, 180 mg of crude, labeled p-fruc- 
tose was dissolved in a few drops of water, and the 
solution was filtered through a microfilter, containing 
some decolorizing carbon, into a 50-ml round-bot- 
tomed flask. The funnel was carefully washed, and 
the combined solution and washings, after the addi- 
tion of a drop of acetic acid, were concentrated under 
reduced pressure to a thick syrup at a temperature 
less than 35°C. The sirup was dissolved in absolute 
ethanol, and the solution was transferred to a heavy- 
walled test tube and reconcentrated under reduced 
pressure. The residue was dissolved in methanol, 
the solution was again evaporated under reduced 
pressure, and the thick sirup was diluted with 5 or 
6 drops of methanol. By cautious addition of a 
methanol-isoamy! alcohol mixture (1:1), and finally 
of isoamy! alcohol alone, the volume was increased to 
approximately | ml, and the solution was brought to 
incipient turbidity. It was then lightly seeded and 
allowed to stand at room temperature. From time 
to time additional isoamyl alcohol was cautiously 
added. At the end of a few days the mother liquor 
was removed, and the crystals were washed in the 
tube with ethanol and dried. To reerystallize, the 
p-fructose was dissolved in methanol, the solution 
was concentrated under reduced pressure, and meth- 
anol and isoamyl alcohol were added as before. 


3.2. Assay of Carbon | of p-Fructose-] 


A 75-mg sample of labeled p-fructose | 
activity of 2.62 ue per millimole was diss: 

ml of water. The solution was _ transf; 
frozen to a portion of the wall of a 200-m! 

flask, previously cooled in a dry ice-aceton 

On a separate portion of the wall, 2 ml of 2 V pos. 
sium hydroxide was frozen. The flask was 

ately placed on a shaking machine,‘ evacua 

filled with oxygen under 10 |b of pressur: 

then shaken for 24 hr at room temperature. | 
procedure, the solutions of alkali and sugar di 
come into contact except in the presence ¢ 
oxygen, and the deleterious effect of alkali upo 
sugar in the absence of oxygen is avoided Aft 
hr, 50 ml of methanol was added and the solyry 
was lightly seeded with dust of crystalline potass 
p-arabonate and allowed to stand for 1 da 
mother liquor was decanted from the crystals, w! 
adhered to the walls of the flask, and the erystals y 
carefully washed with methanol. The 
p-arabonate was recrystallized once in the flask 
dissolving it in a minimum quantity of water , 
adding methanol to the point of incipient turbid 
From time to time additional methanol was adde: 
effect more complete crystallization. After | ¢ 
the mother liquor was removed with a capil 
pipette. The crystals were washed with metha 
and then dissolved in water The solution was 
moved from the flask, filtered through a microf 
containing some decolorizing carbon, and con 
trated to approximately 1 ml in a small test 
under a stream of dry air; methanol was then a 
to the point of incipient turbidity. After erystal 
tion had taken place, the mother liquor was remo 
by means of a capillary pipette, and the product 
washed with 50-percent aqueous methanol. ' 
potassium p-arabonate was recrystallized sev 
times from water, by the addition of methano 
give a product with an activity of 1.32 ye per! 
mole. This value was unchanged by further 
crystallization. Thus the loss in activity accompa 
ing removal of carbon 1 was 49.6 percent. 

It was later found that the potassium p-arabor 
could be purified to constant activity with 2: 
recrystallizations, if the final crystallization 
conducted slowly by allowing an aqueous solu 
in a small test tube to stand undisturbed in a desi 
cator containing anhydrous calcium sulfate and as 
a beaker of methanol. In the course of 1 or 2 ¢ 
large pure crystals formed. The activity of the 
was unchanged by an additional recrystallizatior 


potass 


3.3. Assay of Carbon 6 of D-Fructose-1,6-C" (Carbo 
5 of p-Arabonic Acid) 


For the assay of carbon 6 of the labeled p-fructos 
potassium p-arabonate was prepared and purified ® 
described above, and then treated in the follow 
manner: 


‘A Parr reduction apparatus (shaking device equipped with a 
tank 
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ng quantity of the labeled potassium p- 
e having an activity of 1.18 we per mi'limole 
ted with 1.5 ml of 0.3 N sodium meta- 
and the oxidation mixture, in a glass- 
d tube, was allowed to stand for 3 hr at 
Aqueous sodium bisulfite was 
lropwise in an amount sufficient to cause 
of the iodine that first appeared. The 
was neutralized to the end point of methyl 
by the dropwise addition of aqueous sodium 
nate. To the reaction mixture was added 2.5 
freshly prepared solution containing 35 mg 
done and 35 mg of sodium bicarbonate. 
he solution was carefully neutralized with 
s hydrochloric acid to the end point of methyl 
a voluminous crystallization took place. 
xluet was separated and washed with water. 
ecrvstallization, the material on the filter was 
d in the minimum quantity of methanol; 
ition was passed through the filter and 
with an equal volume of water. After one 
allization, the dimedone compound had an 
of 1.20 we per millimole, which was un- 
d by further reerystallization. Thus, within 
error, the activity of the dimedone 
suund accounted for all of the activity of the 
im p-arabonate. Inasmuch as _ formalde- 
derived, during the periodate oxidation, only 
carbon 5 of the p-arabonic acid (carbon 6 of 
rent p-fructose), 50 percent of the activity of 
eled p-fructose was found at carbon 6. 


mperature 


its of 


3.4. Assay of Carbon 1 of D-Glucose-1-C" 


alkaline oxidation previously deseribed was 
carried out on p-glucose-1-C™ having an activity 
46 we per millimole. In this case the potassium 
rabonate obtained was expected to be inactive. 
salt was recrystallized once in the pressure flask 
lonce ina small test tube by addition of methanol 
solution. Finally, the material in 
ieous solution was filtered, and allowed to crystal- 
slowly, as described in section 3.2. The product 
is obtained showed negligible activity (0.001 ye 
millimole Hence the parent sugar was labeled 
isively at carbon 1. 


al aqueous 


] 


3.5. Preparation of Potassium pD-Arabonate-5-C" 
from d-Mannitol-1-C" 


\ sterile solution consisting of 2 millimoles of 
p-mannitol-1-C™ (having 150 ue of radioactivity 
x of yeast extract, 60 mg of potassium dihydro- 
gen phosphate, and 20 ml of water was incubated at 
(0° C with a freshly prepared culture of A. suboridans. 
The surface to volume ratio was 1.6 cm?/ml.) 

5 hr, 2 ml of a 20-percent solution of zine 
was added, and then an equivalent quantity 
ous barium hydroxide. The resulting precipi- 
is separated by filtration through a funnel 
coated with diatomaceous earth and decolor- 
The filtrate was concentrated to a 


tH} I 


irbon 


volume of about 5 ml, and then frozen on the wall of 
a 200-ml pressure flask The volume with washings 
was about 8 ml Eight milliliters of 2 N potassium 
hvdroxide was frozen on another portion of the wall, 
and the oxidation was conducted as deseribed in 
section 3.2. The resulting solution was diluted with 
5 volumes of ethanol and allowed to stand 24 hr for 
the crystallization of potassium p-arabonate. After 
separation and recrystallization, the product weighed 
214 mg (52.4-percent yield 


3.6. Preparation of D-Arabinose-5-C" 


One millimole (204 mg) of potassium p-arabonate- 
5-C™ was dissolved in 10 ml of water and passed 
through a column containing 10 ml of cation ex- 
change resin; ° this was followed by 20 ml of wash 
water. The effluent was concentrated to a sirup 
under reduced pressure, and the sirup was transferred 
to a reduction tube such as that described in [7]. 
The sugar acid was lactonized by heating it in a 
current of air at 50° C for several days, with occa- 
sional addition of methanol. When the product had 
become entirely crystalline about 1 week), the tube 
was immersed in an ice bath and connected with a 
stirrer, which reached to the bottom of the tube. 
After the addition of 3.2 g of crystalline sodium acid 
oxalate and 10 ml of water, vigorous stirring was 
begun, and 4.6 g of 5-percent sodium amalgam 
pellets was added through the sidearm of the reduc- 
tion tube. Stirring was continued until the amalgam 
was spent, after which the mercury was separated. 
Five volumes of methanol was added; the crystalline 
sodium salts were separated by filtration, washed 
with methanol, and discarded. The alcoholic solu- 
tion was concentrated under reduced pressure to 
about 5 ml, 25 ml of methanol was added, and the 
solution was filtered to remove a second crop of 
nonradioactive sodium salts. The solution was con- 
centrated to remove the methanol, and was neutral- 
ized at ice temperature with sodium hydroxide to the 
end point of phenolphthalein. A small crop of non- 
radioactive salts was removed, and the solution was 
then passed through an ion-exchange column con- 
taining equal parts of cation and anion exchange 
resins. The salt-free and neutral effluent was con- 
centrated to 100 ml under and 
finally was freeze-dried. The residue was dissolved 
in 0.5 ml of methanol, and isopropanol was added to 
the point of incipient turbidity. Crystallization of 
p-arabinose-5-C “ progressed readily, and ultimately 
103 mg of crystalline p-arabinose was separated 
without the addition of carrier (vield=68.6 percent).? 
As potassium p-arabonate obtained in 52.4- 
percent yield from p-mannitol, with the elimination 
of one-half of the carbon 14, the over-all radio- 
chemical yield of p-arabinose-5-C™ was 18 percent 


ee 


Amberlite IR 120, Rohm & Haas Co., Philadelphia, Pa 
* Amberlite IR 100, analytical grade, Rohm & Haas C Philadelphia, Pa., and 
Duolite A-4, Chemical Process Co., Redwood City, Calif - 
’ There has been considerable veriation in yield in this preparation 
chemical yield of 60 to 70 percent may be expected. 


reduced pressure, 


was 


a radio- 
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Measurement of Variations in Atmospheric Refractive 
Index With an Airborne Microwave Refractometer'’ 


Howard E. Bussey and George Birnbaum 


4 microwave refractometer for aircraft 


use is described, 


and some of the sources of 


error in the measurement of variations of atmospheric refrac tive index with this instrument 


are discussed 
wo refractive- 
aspects, resembled changes shown by 
erratically with time 
changing. A rough 


index soundings taken 14% 


analysis of the 


Observations were made up to 10,000 feet on 2 days near Washington, D. C 
hours 
radiosonde 
and place, but was usually greatest where the 
data showed that the 


apart showed changes that, in the coarser 
data. The fluctuation intensity changed 
vertical gradient was 
larger fluctuations occurred over 


distances of several hundred meters, whereas over distances less than 5 meters the fluctuations 


were negligible (the response of the 
tions occurring within 0.5 meter 
cumulus clouds and intense 


Large 


1. Introduction 


rience shows that meter- and centimeter- 

electromagnetic waves are propagated far 

d the horizon [1, la] ? with a field strength much 

than what can be accounted for on the basis 

of standard diffraction theory [2]. To explain this 

discrepancy it has been proposed that the waves are 

red by inhomogeneities in the refractive index 

atmosphere. Consequently, there is need for 

tailed knowledge of the refractive-index structure 

atmosphere, especially with regard to sharp 

ges in the vertical gradient and irregular local 

riations from the mean value. Apart from these 

considerations, observations of refractive-index vari- 

in the upper atmosphe re are of interest in their 

ight in mete orologic al espec ially 

their relation to the water-v vapor 
there 

\lany observations of atmospheric refractive index 
have been made recently with an airborne microwave 
refractometer by Crain and coworkers [4]. <A few 
‘hbservations with an airborne microwave refractom- 
eter of somewhat different design are reported here. 
Except for modifications necessitated by aircraft 
installation, this instrument was essentially that 
leve ‘eloped by Birnbaum [5] and used for observations 

the ground [6, 7]. Special attention was given 

problem of measuring rapid fluctuations as 

is virtually no information available on refrac- 
tive-index inhomogeneities of small extent. 

This paper dese ribes the essential features of the 
mirborne instrument and its performance, exhibits 
elected portions of the data obtained on a flight up 
» 10,000 ft and through clouds, presents the results 
f a crude statistical treatment of the data, and 

isses the experiment in relation to previous 
rological measurements. 


science e, 
because of close 


structure 


Third Radar-Weather Conference 
and at the Joint Meeting of 


his paper were presented at the 
rsity, Montreal, Canada, Sept. 1952, 
IRE, Washington, D. C.., April 1953. 
brackets refer to literature references at the end of this paper 
{ radio-wave scattering by the atmosphere nas been published by 


instrument would have 
increases in 


fluctuations were noted within the 


allowed the detection of fluctua- 
index were observed on entering 
clouds 


2. Description of Equipment 


The frequency change, Af, of a cavity resonator due 
to a refractive-index change, An, of the air contained 
therein is expressed to a very good approximation by 

—Af/f=An. 
As in this work An is of the order of 10~°, it is con- 
venient to introduce a quantity AN defined by 


AN An X10 


where N is called an N-unit. Af is measured|by 
electronically comparing the frequencies of a ven- 
tilated (test cavity and a sealed (reference) cavity. 
These are excited by a frequency -modulated klystron, 
and their outputs after detection vield resonance 
curves that are used to control the width of a rec- 
tangular voltage waveform. This voltage is aver- 
aged, amplified, and recorded. 

The design of the instrument constructed for air- 
craft operation differed somewhat from that de- 
scribed [5, 6] im previous reports. The klystron (type 
2K25) was isolated from the large pressure changes 
encountered in ascents and descents by sealing its 
mount. An available power supply originally 
designed for 60-c/s input voltage operated satisfac- 
torily with the 800-c/s voltage supplied by the 
aircraft. The klystron sweep voltage was synchro- 
nized to the aircraft mains, and consequently, the 
time constants of the electronic circuits were adjusted 
for a pulse repetition rate of 800 ¢/s. 

The sampling cavity * was installed at a location 
believed to be in free air, undisturbed by turbulence 
along the fuselage. It was attached on the end of a 
streamlined boom (fig. 1) 3.7 ft long located 6.7 ft 
from the nose of the aircraft (a Navy type P4Y2). 
The waveguide lead to the cavity was inside the 
boom. The electronic circuits, reference cavity, and 
klystron located nearby were shock-mounted. As 
further protection from aircraft vibration, low micro- 


¢ This and the reference cavity operated in the T Eo, mode at 9,000 Mc 
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10 SECONDS 
750 METERS 


Ficure | 


phonic tubes (1L2AY7) were employed in the input 
stages of the resonance pulse amplifiers Connec- 
tion between the shock-mounted waveguide and the + 
rigidly mounted waveguide in the boom was effected se ii oN 
by a short length of flexible coaxial cable 

The reference cavity was provided with a small 
diameter plunger driven by a micrometer whose 
scale was calibrated in terms of cavity frequency 
change. By means of this plunger the refractometer 

The response time of the present eq 


was calibrated, also the recording pen was reset reproduced. 
ment, however, was found to be limited by tha 


the recording milliammeter.© The response of 
recorder was essentially flat out to 20 e/s, a 
decreased by a factor of one-half at 45 ¢/s. 


Ficure 2 Two segments of the refractive-index reco 
/ } 


while in level fight in very homogeneous air at abo 


whenever it approached the off-scale position. The 
latter adjustment permitted the use of a sensitive 
scale to determine fine detail, vet preserved the much 
larger changes in refractive index that occurred in a 
sounding While making atmospheric observations a const 
chart speed of 2.5 mm/sec was used. Consequer 


Invar of high purity that had been subjected to a 
oscillations more rapid than 10 ¢/s would not 


special heat treatment [8] was used in the construc- 
tion of the test cavity in order to minimize thermal individually resolved. They would, how 
expansion. The treated metal had «a measured broaden the trace if they were continuously pres 
coefficient of linear expansion of 0.3 10~ + as the recorder could respond to higher frequen 
and the finished cavity had the same coefficient The trace width was normally '; mm, which 

because the change in the refractometer reading due equivalent to 0.05N on the most sensitive ra 
to heating of the cavity was 0.3.N/° C However, used, 7.N full scale. On the basis that an 80-p: 

after installing the cavity in the aircraft the coeffi- broadening of the trace was easily visible and t: 

cient increased to 0.510 C. It is felt that account of the decrease in recorder respons 
this was due to stressing by the aluminum supports. higher frequencies, it was estimated that the n 
The reference cavity available from previous work mum amplitude of index oscillation necessary 
Was temperature compensated to within0.5*10~°/° CC detection at various frequencies was as indi 
by using brass end plates suitably recessed into an in table | 

untreated invar (a=3 10 C) barrel rasie 1. 

Possible changes in air density arising from the 
rapid passage of air through the cavity were investi- 
gated at speeds up to 170 mph. The test cavity 
was mounted in the NBS wind tunnel and changes 
of air density in the cavity were deduced from meas- 
urements with the refractometer. Within the exper- 
imental error of 0.5 to LN caused by constant 
addition of fresh air to the tunnel, the densities of 
air within tunnel and cavity were alike 

The smoothing of air-temperature fluctuations due The length there is the distance traversed b 
to contact of air with the cavity walls was determined aircraft (whose speed was 75 m/sec) during a 
in previous laboratory tests for air speeds up to 25 cycle of the oscillation. It indicates the app 
mph. At this speed the amplitude of air-tempera- Mate “size’’ of the index inhomogeneity that » 
ture fluctuations was reduced by 20 percent; but at gtve rise to the frequency in question. 
aircraft speeds the reduction should be negligible. The smallest refractive-index fluctuation that 

The ability of the refractometer to record rapid be reliably observed is limited by noise origina 
changes in index clearly depends on the repetition within the equipment. In the laboratory a 1 
rate of the klystron sweep voltage. Thus for a ™um noise amplitude of 0.05N was observed 
sweep rate of } equal to 800 ¢/s in the present case), this occurred in intervals of the order of 1 to 2s 
oscillations at a frequency of f/2 should be detected There were virtually no current variations visi) 
and the shape of those at //6 or less should be well 


By strongly modulating the screen voltage of one of the amplific 
the refractometer rectangular pulse, and hence output current 


vary sinusoidally 
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rvals or shorter. In flight it was not cavities as the temperature environment changes 
t to seal the ventilated cavity in order to The observed temperature change in the ascent to 
record of equipment noise. Instead, an 10,000 ft was large enough to give an indicated 
of this noise was estimated by assuming change of several N-units. Some variation might 
one was responsible for the variations in occur also from drift in klystron center frequeney 
hest section of the data. Figure 2 shows with temperature, as the klystron was not thermo- 
records obtained in very homogeneous air statically controlled. In addition, pressure charges 
ft. The deviation from the mean in \-see during ascent might shift the frequency of the sealed 
was about 0.02N and in l-see intervals cavity slightly Although these effeets could cause 
rp errors of perhaps 5 percent in the sounding, they 
rt was made to determine the longtime do not influence the reliability of the observations 
f the equipment. This depends in large — of fluctuations. 
however, on the thermal expansion of the 
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3. Results 


The observations described in this section were 
obtained over the Chesapeake Bay within 12 miles 
of North Beach, Md., on August 19, 1952, from 1300 
to 1600 EST. The weather was warm and humid. 
There was a 15-mph north wind up to about 8,000 
ft, and above this level there was a 25-mph west 
wind. The radiosonde temperature profile indicated 
a dry adiabatic lapse rate (neutral equilibrium) up 
to 4,000 ft and a moist adiabatic lapse rate above 
4,000 ft. The presence of scattered cumulus clouds 
with bases at 4,000 ft and tops at 12,000 to 14,000 ft 
was evidence of large-scale vertical convection 
processes 


3.1. Description of Data 


The horizontal speed of the aircraft was kept to 
about 75 m/sec (170 mph) and during a sounding 
the ascent or descent rate was held to about 2. 
m/sec. Thus the slope of the inclined path was about 
one-thirtieth. In_ figure are shown sections of 
strip charts, each of which is 126 see in duration. 
During this time the aircraft traveled 10 km hori- 
zontally, and, if making a sounding, changed altitude 
by 300 m 

Although the atmosphere at 10,000 ft was for the 
most part smooth and homogeneous, as shown in 
figure 2, there was a short time when it was not, as is 
illustrated in figure a. This figure shows the 
occurrence of slow variations of appreciable ampli- 
tude, but with virtually none of the rapid fluctuations 
that were usually observed at lower levels. 

Figure 3, b, illustrates the pronounced and wide- 
spread uniformity of refractive index that was 
encountered from 8,500 to 8,300 ft (during the 
descent), and which persisted for the most part up to 
10,000 ft the highest altitude explored. 

Figure 3, c, shows a very large increase in index, 
namely, 25.V, between 8,300 and 8,100 ft. (The pen 
went off scale during this increase and was reset at 
the place marked RS.) This region of fast transition 
was caused by a widespread and persistent stratifica- 
tion of the atmosphere. It was observed during four 
ascents or descents within a half-hour period and was 
always encountered at the same height to within 
+100 ft. This apparent boundary between two air 
masses was visible as a sharply defined haze line, 
with comparatively clear air above. This would 
indicate, in agreement with the large variation of 
index, that the air was moist below the transition 
and dry above it. 

Another striking feature of this record is the 
smooth and very rapid increase of 10N at 7,900 ft. 
This increase, which represents one of the sharpest 
boundaries encountered, occurred within 0.7 sec, 
during which time the aircraft traveled 50 m hori- 
zontally and descended 2m. The index then slowly 
returned to about its former value in a horizontal 
distance of 2,500 m. 

At 7,000 ft the plane was leveled for horizontal 
flight. Fluctuations of various durations and ampli- 
tudes occurred (fig. 3, d), but most interesting of all 
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was a well-defined refractive-index variation 
extended about 5 miles and had a mean amp 
about 12N. 

Figure 3, e, was 
7,000 to 6,000 ft. 

An example of a record in which the full 
sensitivity was changed from 33 to 7N is show 
figure 3, f, although in general the index varia; 


obtained in descending 


», 
were too great to permit the use of the more sensi: 
scale. No perceptible broadening of the stylus 
occurred on changing sensitivity, indicating 
continuous fluctuations with durations of 0.1 to 
see (7.5 to 1.5 m) were of amplitude less than 0) 
The smallest fluctuations that are plainly vis 
have amplitudes of 0.1 to 0.2N and durations of 
sec. There are many fluctuations of 1 to 2\ 
durations of a few seconds. 

At 3,000 ft the fluctuations (fig. 3, g) were defin 
of greater amplitude than at the higher elevat 
The strikingly large-amplitude and long-dura: 
fluctuations near the end of the record are indica 
of large and rather homogeneous eddies of air w! 
refractive index differs markedly from the n 
As will be indicated later, these variations as wi 
those of smaller amplitude arise chiefly from 
tions in vapor density. 


9 
»,. 


3.2. Refractive Index of Clouds 


Refractive-index observations were made at 7 
ft while flving through some fair-weather cun 
clouds. These had bases at roughly 4,000 ft 
tops in the neighborhood of 12,000 to 14,00 
The air temperature inside and outside the « 
was approximately the same, 13.5° C. On ent 
a rather well-developed cloud the refractive 
increased by about 40N and then fluctuated 
lently (fig. 4, a). The increase that can be se: 
the figure amounts to 16.V, but while the trac 
off scale (marked A) the frequency of the refer 
cavity was changed by an amount equivalent to 
in an attempt to set the pen on seale. Ther 
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ddenly decreased and went off scale as the 
nse part of the cloud was passed The trace 
tat points Band C. By taking into account 
lacements of the tuning plunger, it was 
t after leaving the cloud the refractive index 
to its original value. Changes of 6 to 10.N 
distance of about 25 m were observed in the 
Further evidence of the great turbulence 
sted there were the bumps experienced by 
aft. By way of contrast, it may be men- 
that a moderate bump was encountered in 
with no noticeable change in index. 
itering another cloud not as well developed 
rst, the index increased by 35.N; and on leav- 
reased by this amount (fig. 4, b In this 
full-scale sensitivity was 135.N. Inside the 


everal changes as large as 20.NV were noticed. 


irkably sharp refractive-index gradients were 
d in both clouds; changes of 12 to 15N 


a distance of 12 m occurred in three or four 


s and of 6N within 2 m in one instance. 
circling at 10,000 ft, where the refractive 
vas quite uniform, it was noticed that a 


A smoothed ascent sounding, and a descent sounding 
taken about 1‘4 hours later. 


e descent curve occur due to interruption of the descent for hori- 
rations. The adjectives describing the rapid fluctuations (appli- 
lescent curve only) represent amplitudes approximately as fol- 
ry small, 0.25N or less; 8, small, 0.5N; L, large, 1.0N; VL, very 
more 


definite increase occurred as a nearby cloud was 
passed. The distance of closest approach to the 
cloud was 1 or 2 miles 


3.3. Refractive-Index Profiles 


To obtain the variation of refractive index with 
height, the original data were smoothed by obtaining 
visually sucessive 10-see averages. Curves plotted 
from such values are shown in figure 5. Included 
there are rough estimates of the mean amplitude of 
those fluctuations less than 5 see in duration. On 
comparing these estimated amplitudes with the 
gradient, one notices the not surprising result that 
regions of change of gradient are associated with 
large fluctuation amplitudes 

One of the curves in figure 5 is based on data taken 
in ascent and the other on data taken in descent 
about ls hrlater. During this interval the refractive 
index clearly increased at all levels above 5,000 ft. 
For comparison, refractive-index profiles computed 
from radiosonde data are shown in figure 6. These 


data’ were obtained in the vicinity of Silver Hill, 


Md., in the morning and evening of the day of the 
flight. The increase in index during this longer 
interval is somewhat greater than that noted above 
This increase as shown by the radiosonde data must 
have been due almost entirely to changes in water- 
vapor density, as the temperature and pressure 
profiles scarcely changed. 
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August 19, 1952, calculated from radiosonde ascents 
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4. Data Analysis 


presents the results of a rough 
statistical treatment of the data upon the 
following considerations. The atmosphere may be 
regarded as an assembly of air parcels whose refrac- 
tive-index values differ from another and in 
general from the mean value. Within any large 
parcel there may be smaller parcels of smaller index 
variations and soon. The refractive-index amplitude 
of a parcel was taken as the peak positive or negative 
deviation from the local mean and the extent or size, 
in the direction of flight, as the interval between 
consecutive intersections of the record with the local 
mean. A rough estimate was made of the percentage 
of time that parcels of various sizes appeared to be 
present. The duration of the local running mean 
was adjusted in proportion to the size of the parcels 
to be analyzed. The results are expected to be 
roughly related to Fourier amplitudes and periods as 
might be obtained by usual statistical methods. 

The results of the data analysis are summarized in 
table 2, which shows that both the amplitude and 
percentage of time present increase with the size of 
the parcel. In a few cases, where level flights in a 
given direction were made, sizes greater than the 
6,000 m indicated in table 2 were observed Indeed, 
the presence of air parcels some 10 to 30 miles in 
extent has been detected by aircraft measurements 
of temperature and vapor pressure |9, 10] 
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The 1,000-m range of sizes at altitudes above 
8,500 ft were not included in table 2. Variations in 
the altitude of the aircraft over such distances could 


have caused the small refractive-index va 
observed. 

For parcel sizes of the order of several | 
meters, the variation of amplitude with hx 
shown in figure 7. In this particular case 
certainly no clear systematic variation. It sh 
noted that the amplitudes for the region | 
8,000 ft in table 2 are averaged over the wi 
ations indicated by figure 7. 

An important finding of this analysis was | 
that the amplitude of the refractive-index 
geneities of small extent was extremely small 
the absence of broadening of the trace whi 
sensitivity was 7.V full seale, it was conelud: 
sizes in the range of 0.5 to 5 m, if present, had 
tudes no greater than 0.04N. There wer 
10 instances in the flight where this more se 
scale was used. In addition there was an initial | 
flight (made on August 15), with most of th 
recorded at 10,000 ft, which also showed the al 
of trace broadening within the limits specified abo 
As for the remainder of the present flight, taken \ 
a full-scale sensitivity of 33.N, from the absence 
trace broadening it can only be stated that 5-m s 
if present, had amplitudes less than 0.1.N at all tir 


5. Discussion 


5.1. Causes of Refractive-Iadex Fluctuations 


For the mean conditions of temperature, va; 
pressure, and total pressure encountered in the flig! 
the change in refractive index, An, is approximat 
related to the change in temperature, A7(° C 
to the change in vapor pressure, Ae(mm Hg) by’ 

An X10°=AN=6de— AT. 
Rapid pressure fluctuations are too small to influe 
AN. Radiosonde records * taken at Silver Hill, M 
on the day of the flight were examined, and it ) 
found that AJ fluctuations for sizes that the ra 
sonde can best detect—50 to 100 m—were 
less than 0.2 deg C, thus accounting for only a fraet 
of AN for these sizes (see table 2). It 1s proba! 
that for other sizes also the changes of index w 
mainly due to changes of vapor pressure. This \ 
is supported by the measurements of reference [4 


usua 


5.2. Clouds 


The large increases in N observed in the clouds 
reasonably accounted for on the basis of changes 
water-vapor content. The uprising saturated au 
the clouds at a temperature of 13.5 ° C should hay 
refractive-index value about 70.N greater than 1! 
were dry. As shown by radiosonde data, the rele 
tive humidity outside the cloud at 7,000 ft was alx 
40 percent. The observed increases of 35 to 40.\ 

® This result was obt 
e/T?, using representative values 
mm He 


* Provided by the U.S 
allow parcels of dimension 


ned from the expression 1)108— 103.5 
of the variables P is the t 


Wether Bureau 
in the vertical direction 


The time constant (4 
of Sto 100m to 


176 





“d abo 
ken \ 
bsene 


-M s 


all times 


itions 


e, Val 


he flig! 


( 


by : 


influer 


fill, M 


id it was 


he rar 
* Usua 


i fract 


probal 


dex W 
This \ 


ice 4 


louds 
anges 
ed all 
id hav 
han i 


the reli 


as ab 
Oo 10.\ 


Ximat 


ds are therefore reasonable. The liquid drop- 
hese clouds would contribute a small amount 
idex—less than LV. An artificial effect due 
sition of liquid drops on the cavity walls does 
ear likely because of the rapidity of the ob- 
ariations. The existence of large amplitude 
ons within the cumulus cloud may be evi- 
the process called entrainment [11], in which 
ding drier air mixes with the uprising satu- 
r of the cloud. 
ndex inhomogeneity in the clouds was so 
hat it was wondered whether appreciable 
turn could result. The back-seattering cross 


per unit volume is given by An?/Sal [15] 
the seale of turbulence, is large compared to 


r wavelength. From the data it was esti- 
hat 7 was of the order 100 m, and An? was 
These values vield a cross section of 410 

For 10-cm-wavelength radiation, this re- 
bout 100 times larger than that expected due 
leigh seattering from rainless clouds in 
he mean drop diameter is 0.001 em, but even 
mounts to only a very small scattering cross 
capable of detection only at a short range 
der of 1 mile. 


5.3. Soundings 


iscent and descent soundings (fig. 5) differ not 

as to the aspect of the coarse cradient but also 
the amount of fine structure The ascent was 
mmediately after low-lying stratus clouds 
dissipated by the sun’s heating and may have 
el place before the onset of strorig vertical 


The subsequent descent shows greater 


homogeneity, such as might have been caused by 


espread vertical mixing 
It should be realized that the various changes of 
ulient in a sounding may be due to three causes 
a rather permanent stratification,'’ (2) a tempo- 
v fluctuating gradient, and (3) an artificial 
t that occurs because the angle of ascent is small, 


hus allowing horizontal gradients to give rise to 


rent vertical gradients. The existence of per- 
lavers can, of course, be confirmed by re- 

ed soundings as was done in connection with the 
boundary at 8,300 ft. The large gradient 
ved in ascent at 3,500 ft may also have rep- 
ited a laver boundary, but this was not checked 


5.4. Index Parcels of Small Size 


present work it was found that the fluctua- 
nplitude associated with an index parcel 
d with the size of the parcel. In fact, the 
des of 5-m parcels (and smaller) were found 
egligible. A similar finding was reported by 

al. [4], although their observations did not 
to as small a parcel size as observed here. 


2] may have det 


These results, although of limited scope, lend some 
support to Megaw’s [14] assumption, made in con- 
nection with his theory of atmospheric scattering of 
radio waves, of a (Kolmogoroff spectrum of parcel 
sizes terminated at 10 em by viscosity on the small- 
size end. Booker and Gordon [15], in dealing with 
radio-wave scattering, used a scale or mean parcel 
size of 10 em with an amplitude of 1.N, but no such 
parcels were observed in the present experiments 


6. Conclusions 


Although a limited amount of data was taken, the 
results served to demonstrate the operation of the 
instrument and to give some notions about the 
refractive-index structure of the atmosphere. It 
was found that the character of the fluctuations 
changed erratically with position and time. <A 
persistent feature of the two flights was the absence 
of very rapid fluctuations. Quantitatively expressed, 
it was found that fluctuations that occurred in 
distances of from 0.5 to 5 m had an amplitude less 
than 0.04 107° in n (or O.O4N The larger ampli- 
tude fluctuations occurred over distances of the 
order of 1,000 m, which interestingly enough is in 
general agreement with the mean seale of turbulence 
of 300 m for the atmosphere obtained by Megaw 
14] from stellar scintillation data, and with hori- 
zontal dimensions of 2,000 to 3.000 m of some of the 
parcels observed by Bunker [16] 


The measurements described here were made in 
cooperation with the Naval Research Laboratory 
and the aircraft was furnished by that laboratory 
We express our gratitude for the cooperation extended 
to us by the members of the Wave Propagation 
Research Branch, and especially to R. BE. Katz, who 
aided in several phases of the experiment and to 
H. D. Metke, who piloted the aircraft. A. Hanyvok 
constructed the electronic equipment. We also 
thank H. Lyons, K. A. Norton, and J. W. Herbstreit 


for their interest 
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Studies in the System Magnesia-Silica-Water at Elevated 
Temperatures and Pressures 


Elmer T. Carlson, Richard B. Peppler, and Lansing S. Wells 


The reactions between magnesia and silica under hydrothermal and pneumatolytic 


conditions were investigated. 
of the reacting oxides. 
microscope to have a fibrous structur>. 


Chrysotile and tale were produced, depending on the ratio 
The chrysotile was submicroscopic but was shown by the electron 
Chrysotile and tale were likewise 
mixtures of magnesium carbonate and silica. 
action of silica in the vapor phase on solid magnesia. 


prepared from 
The same products were obtained by the 
A detailed study of this method was 


made, and the effects of temperature, pressure, time, pH, and oxide ratio on the amount 


of silica transported were determined 
with temperature, pressure, and pH 


1. Introduction 


The system MgO-SiO,.-H,O has been the subject 
extensive study by a considerable number of 
estigators over a period of many years. This 
videspread interest is in part due to a desire to 
explain the mechanism of the formation in nature of 
the various hydrated magnesium silicate minerals, 
particularly serpentine. Probably a stronger motive 
for much of the research has been the hope of 
developing a method for the practical synthesis 
of serpentine in its fibrous form, chrysotile, which is 
most important of the asbestos minerals. 

The various terms applied to the serpentine 
minerals are not used in the same way by all investi- 
In this report no attempt will be made to 
justify the terminology employed, but the mineral 
names will be used in accordance with what appears 
common usage. The word serpentine 
will be applied in general to the mineral or to 
synthetic preparations having the approximate 
composition 3MgO-2SiO,-2H,O. The terms chry- 

« and antigorite will be applied respectively to 
the fibrous and the platy crystalline forms of ser- 
pentine. These two forms presumably differ in 
their internal structure, but the X-ray patterns that 
have been published have many features in common, 
and often are not free from ambiguity. 

Magnesia and silica have been found to combine 
readily in the presence of water vapor, especially at 
elevated temperatures and pressures. Among the 
products obtained is a compound having the chemical 
composition of chrysotile, but never its fibrous habit, 
far as can be seen with the light microscope. 
Bowen and Tuttle [1], in their study of the system 
Mg0-SiO.-H,O, showed by differential thermal 
and other means, that this actually is 
ile. Electron photomicrographs showed that 
ystals were definitely rod shaped. These 
s also showed that chrysotile may be formed at 
atures below 500° C, but not above, and at 
pressures from 2,000 to 40,000 Ib/in?. 
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The overall rate of transport was found to increase 


One of the ultimate objects of the present work was 
the preparation of chrysotile fibers of larger size. 
This attempt proved unsuccessful, but some of the 
information obtained appears of sufficient value to 
be recorded. 

As a preliminary step in this investigation, an 
exhaustive search was made of the literature per- 
taining to asbestos, as well as to the system MgO- 
SiO.-H.O and related subjects. A bibliography of 
selected references has been compiled.” 


2. Experimental Studies 
2.1. Solubility of Magnesia 


As a starting point it was considered desirable to 
know the solubility of the reactants at elevated 
temperatures. The solubility of silica in water has 
already been investigated thoroughly over the range 
from 0° to 425° C, and at pressures up to 1,250 atm. 
The available data being adequate, no further work 
was done in this field. 

The situation with respect to magnesia is quite 
different. A considerable number of investigations 
of the solubility of magnesia at and near room tem- 
perature have been made, but the data vary over a 
wide range. It has peen shown that the solubility of 
magnesia varies with the temperature at which it was 
prepared, and varies also if it was prepared at a 
given temperature from different materials. Gjald- 
baek [2], in 1925, concluded that there is an unstable, 
soluble form of Mg(OH). which first forms and is 
quickly converted to a stable, less soluble form 

There appears to have been only one investigation 
of the solubility of Mg(OH), in the absence of CO, 
and at temperatures above the melting point of 
water. This is the work of Travers and Nouvel [3] 
in 1929, wherein the solubility of Mg(OH), was 
investigated in the range from 0° to 200°.C. The 
description of their method is very brief, and nothing 
is said about filtration of high-temperature solutions. 
It therefore seemed advisable to investigate the 


*A limited number of single ~ pies are available from the authors, 
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was an initial supersaturation lasting about 
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if magnesia in 


solubility of magnesia over as wide a temperature 
range as possible 

The solubility 
same MgO that was used in subsequent synthesis 
experiments. It was of laboratory quality and had 
been previously ignited for 1 hr at 1,000° C. It 
was found to contain 0.2 percent of extraneous 
soluble salts. A series of solubility determinations 
were also made on another portion of this same 
MeO that had been extracted with water and that 
had not been previously ignited. The values ob- 
tained with this MgO were slightly but consistently 
lower than with the ignited and 
unwashed MgO 

The water used in each case was ordinary distilled 
water, the CO, content of which had been decreased 
by boiling and bubbling CO,-free air through it. A 
few determinations were made, using water that had 
not been treated in this fashion, and the solubilities 
thus obtained were consistently slightly higher. 

Solubility determinations were made at seven 
different temperatures in the range 30° to 250° C. 
The determinations at temperatures under 100° C 
were made in rubber-stoppered polyethylene bottles, 
that were stored, with occasional shaking, in elec- 
trically heated ovens maintained at the desired 
temperature. Portions of these solutions were taken 
from time to time and analyzed until a value was 
obtained that did not change for several days. The 
MgO was determined by titration against standard 
HCl, using methyl red as an indicator 

For the determinations at temperatures above 
100°, use was made of an autoclave of special design 
described in an earlier paper [4]. This autoclave 
was of stainless steel, with a capacity of 2 liters, and 
was fitted with a Munroe filter, which permitted 
filtration of the solution before the autoclave was 
opened. Aliquots were analyzed as described above. 

The data are presented in table 1 and plotted in 
figure 1. The solubility was found to decrease 
steadily from 2.1X10~* mole per liter at 30° to 
1X10-* mole per liter at 250° 


determinations were made on the 


those obtained 


continuously on a suitable mechanical tumbler, ; 

duration of the period of supersaturation was grea 

reduced but the magnitude remained about the san 

Because the actual shape of the solubility-versys 
time graphs is a function of technique used, they ; 

not reproduced in this paper. 

At the lowest temperatures there was an appr 
able difference between the data obtained for ¢) 
washed and for the unwashed ignited magnesia 
the temperature increased, the difference decreas 
accordingly, at temperatures above 100° C, only o 
kind, the ignited MgO, was used. 

In the temperature range 100° to 250° C 
solubility of magnesia steadily decreases (fig 
The smallest value, 1 10~° mole/liter, is the limit 
usefulness of this method of analysis, for thereafi 
the area of confusion in the end point is of the sa 
order of magnitude as the titration 

The data of Travers and Nouvel [3] are als 
plotted in figure 1. It may be observed that th 
data are of the same order of magnitude as thos 
obtained in this investigation, but consistent 
lower. Because these investigators do not deser 
the magnesia they used nor what precautions, if am 
were taken to remove CO, from the solutions, o1 
cannot comment on this difference 

For the purpose of this investigation it was n 
considered desirable to extend these solubility det 
minations to higher temperatures. It is appare! 
that the solubility of Mg(OH), will become immeas 
urably small before even the minimum temperatu’ 
of the subsequent synthesis experiments is reach 


2.2. Synthesis Experiments 
a, Materials and Methods 


For most of the experimental work the magnes 
used was reagent quality magnesium oxide. This 
was found to be hydrated to a considerable, # 
variable, extent. Calcium in soluble form we 
found to be present to the extent of 0.17 percent 
expressed as CaO. No attempt was made to deter 
mine the actual calcium compounds present 
was not considered necessary to purify the magnes# 
further by extraction with water. Failure to dos 
however, gave rise to some confusing results, to 
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| later. The magnesia was used either in 
tially hydrated state or after ignition at 
temperatures. 
few experiments, reagent quality magnesium 
te was used. Analysis indicated that it was 
the basic carbonate, hydromagnesite 
).. Me (OH )o-3H,O) 
i] forms of silica, including silica gel, pow- 
nt, and quartz, were used. For the majority 
<periments the silica was present as reagent 
cic acid. Ordinary distilled water was used. 
sure bombs of the Morey tvpe served as reac- 
These were made of stainless steel, and 
iverage capacity of about 18 ml. They were 
n vertical-tube resistance furnaces with pyro- 
control, using iron-constantan couples. In 
the fluctuation in temperature over a 
of a week was within a 4-deg range. Pressure 
ther independently controlled nor directly 
d. The values given for pressure in tables 
d 4 were calculated from the temperature and 
o of the quantity of water to the volume of 
ssel. The charge was placed in a small plati- 
icible supported above the bottom of the 
For the vapor-transport experiments, only 
wcnesia was placed in the crucible, the silica 
put in the bomb cavity beneath. In these 
he crucible was left uncovered, and as a rule 
was added directly to the magnesia. Ex- 
here otherwise stated, the quantity of water in 
mb was approximately the amount that would 
ompletely converted to steam just below the 
cal temperature. 
removal from the furnace, the. bomb was 
n air, in front of a fan, and opened. The 
charge was then dried at about 110° C, and examined 
under the polarizing microscope. X-ray diffraction 
analysis and other tests were then made, if desired. 
\-ray patterns were made on an X-ray Geiger- 
ounter diffractometer, using copper K, radiation. 


sels 


Experiments with Mixtures of Magnesia and Silica 


Several series of experiments were conducted on 
timate mixtures of magnesia and silica, both in the 
2 molar ratio required for serpentine and in the 
tratioof tale. Data for a few of these experiments 
presented in table 2, the rest having been omitted 
they contributed no additional information 

i significance 
The reaction products always appeared amorphous 
when viewed under the microscope, except that there 
was often some birefringence, especially in products 
rmed at the higher temperatures and in the more 
mixtures. The index of refraction was 
1.56 in all cases. By means of X-ray 
tion patterns it was shown that crystalline 
il was produced. After 1 day at 200° C, a 
of magnesia and silicic acid in 3:2 molar 
vave a very weak diffraction pattern with 
bands. Higher temperatures and longer 
periods resulted in stronger patterns. After 
at 250°, a strong pattern was obtained, and 
results were found at higher temperatures. 
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riments on mixtures 
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te ind by analy 
Or Me(OH Hy, 


In these cases the pattern was essentially identical 
with one given by a specimen of natural chrysotile 
asbestos 

Similar results were obtained when magnesium 
carbonate was substituted for magnesium oxide in 
the mixtures. The silica displaced the CO, from the 
carbonate 

Mixtures of magnesia and silica in 3:4 ratio 
produced tale at all temperatures tried 

In a series of experiments with the 3:2 mixture, 
very small amounts of various other substances were 
added. These included compounds of various metals 
that are known to occur in natural serpentine, also 
ammonia, ammonium volatile acids, and 
NaOH. Inno case was there any evidence of fibrous 
structure in the product. 


salts, 


c. Vapor Transport Experiments 


It has been shown that silica is appreciably soluble 
in water-vapor, both above and below the critical 
temperature. Synthesis of serpentine through the 
action of vapor-borne silica on Mg(OH), has been 
reported by Syromyatnikov [5] and van Nieuwenburg 
[6]. It appeared desirable to make a study of this 
method of synthesis, with special attention to the 
effect of temperature and various other conditions 
on the rate of transport of silica. 

For these experiments a weighed quantity of 
MgO was placed in a small open crucible supported 
near the top of the bomb chamber, and the silica 
and water were placed in the bottom of the chamber. 
After the desired period of heating, the crucible and 
contents were removed, dried, ignited at 1,100° C 
and weighed. The increase in weight was ascribed 
to SiO.. In one case the silica was determined by 
actual analvsis, with good agreement as to results 

An attempt was made to evaluate separately the 
effects of relative proportion of the oxides, tempera- 
ture, pressure, and pH on the rate of transport of 
silica to the magnesia. In order to accomplish this, 
it was first necessary to determine whether the 
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results were reproducible. This was found to be 

the case, within reasonable limits, provided all 

‘xperimental conditions were fixed. Results ob- 

tained in one furnace did not always agree with 

those obtained in a different furnace. Specifically, ‘isotherm. Pressure 274 atm 

two of the furnaces were found to give results differing : ; 

widely from all the others. This is believed to be - 1 ames 

due to a difference in temperature gradient In the ' 0. 12 08S; 

majority of the furnaces the temperature at the base 978 a 0.1 asst 

of the bomb was found to be about 2 deg higher ; ; ~ : + 

than that at the top. The two exceptions mentioned 113 

above appeared to have a slight gradient in the 

opposite direction. In these, the transport of silica 

was negligible. The data recorded below represent 

experiments in five furnaces giving reasonably 2 0.7 0. cas 
(30 1001 


concordant results. ves 1385 
1037 f 1553 
oa7To 31 1583 


d. Effect of proportion of oxides. : 


It was further observed that, within the limits the. 


A temperature of 380° C was arbitrarily chosen these experiments, the identities of the reactio w 
and a series of determinations made, using the products of the magnesia and the transported sili ond 
apparatus and techniques described earlier. Using were not a function of the initial proportions 
three initial proportions of silica to magnesia, namely,  gilica and magnesia. These reaction products 
10 to 1, 5 to 1, and 1 to 1, and in all cases, 4 g of discussed later. 
water in the bombs, the amount of silica transported Effects of time and te m perature Four temperatures 
to the magnesia was determined at various times. were arbitrarily chosen, namely, 360°, 380°, 40 
The results obtained are shown in table 3 and and 420° C, and determinations were made of t! 
figure 2. silica transported to magnesia after 3, 7, 14, and 2 

The spread of the values shown in figure 2 is days at each of these temperatures. 
within the limits of experimental error, and it is The initial proportion of silica to magnesia was 
therefore concluded that, providing a certain mini- 10 to 1 in each case, and the bomb contained 4 ¢' 
mum amount of silica is present, the transport of water. Because the amount of water used was co! 
silica to magnesia is independent of the relative stant, the pressures that obtained in the bom!» 
amounts of silica and magnesia. The amount of varied with the temperature and also in the case | 
silica used in all cases was suffleient to cover the the posteritical temperatures with the ratio of t! 
bottom of the bomb. It may be that the available volume of water to that of the bomb cavity. T! 
surface area of the silica is a significant factor, but bombs used were fairly uniform in capacity, av 
this was not investigated. It was observed that the for this work were considered identical. Each ise 
smaller the absolute amount of silica initially present, therm then had its own characteristic pressul 
the greater was the weight percent of it that was These pressures were, for the four temperatures | 
transported under the same conditions of temper- 360°, 380°, 400°, and 420°, respectively, 189, 28° 
ature and time. This is to be expected 274, and 312 atm. 
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Silica transport at va s temperature 


data obtained from these experiments are 
ted in table 3 and figure 3. In figure 3 it may 
observed that the amount of silica transported 
ses With increasing time at each temperature 
that it increases with increasing temperature 
avy further be observed that the three post- 
isotherms are of the same type, showing 


rapid early transport of silica relative to the trans- 


at the one precritical temperature. The iden- 
the reaction products between theemagnesia 
transported silica is not a function of temp- 

in this range. These will be discussed later. 
Efect of pressure. A temperature of 420° C and 
time of 14 days were arbitrarily chosen to study 
(feet of pressure on silica transport. The same 


apparatus and techniques previously described were 
sed. The pressure was varied from 15 to 600 atm 


varving the degree of filling of the bomb with 
The actual values of the pressure under 
varving circumstances were obtained by inter- 
tion in the tables given by Van Nieuwenburg 
Blumendal {7}. It was found that at pressures 
30 atm, in the particular equipment that was 
the transport of extraneous material from the 
and metal gaskets made it impossible to 
quantitative data. For this reason, data 


‘btained at pressures above 330 atm are not included 


this report 
results of these determinations are presented 
ble 4 and figure 4. Referring to figure 4, it 
observed that the amount of silica trans- 
at 420° C increases with increasing pressure, 
ly at first, then very rapidly in the pressure 
200 to 300 atm, after which it begins to 


distribution and proportions of the reaction 
ts depend on the amount of silica transported 
sequently on pressure, temperature, and time, 
dentity of these products is not a function 
variables in the ranges investigated. 
reaction products were examined with the 
iphie microscope, and X-ras powder patterns 
de of samples taken from different depths in 
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Ficure 4 ure on silica transport at a constant 


tem perature of 


Effect of press 
120° C and a contsant time of 14 days 


TARLE 4 Effect on silica transport of varying the pressure 
at consiant temperature of 420° ¢ and constant time of 14 


days 


Experi Weight 


ment of silica 


Only 


the platinum thimble containing the charge. 
three compounds were formed, namely, chrysotile, 
tale, and brucite 

In general, tale predominated in the silica-rich 
surface of the reaction product, and its proportion 


increased with increasing amount of silica trans- 
ported. The chry sotile occurred as microcry stals 
visible under the electron microscope only, but 
giving a sharp X-ray diffraction pattern. Prucite 
always predominated in the bottom half of the re- 
action product and sometimes also in the upper 
half, depending on the amount of silica that had been 
transported. In the investigation of effect of pres- 
sure at 420° C, it was found that the proportion of 
tale and chrysotile in the bottom half of the reaction 
product was greater at lower pressures. 

Effect of pH. The effect of increasing the pH on 
the transport of silica to magnesia was investigated 
at 420° C and a constant period of 14 days. In each 
case the proportion of silica to magnesia initially 
present was 10 to 1, the absolute amounts being 1 
and 0.1 g, respectively, and the amount of water 
present was 4 g. 
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The pH was not measured as such > was varied 
by adding increasing proportions of NaOH to the 
silica initially present. The experiments were other- 


wise identical with those previously described. It 


not 


was found that the transport of silica was 
silica 


affected until the molar ratio of NaOH to 
initially present was about 0.15 and that the amount 
of silica transported thereafter increased with in- 
creasing NaOH. When the initial NaOH to silica 
molar ratio was 0.32, the amount of silica transported 
had increased about 64 percent over what was trans- 
ported under the same conditions without any 
NaOH being present 

The NaOH appeared to have no effect on the re- 
action product other than what could be expected by 
increasing the amount of silica transported, that is, 
increasing the proportion of tale 
of experimental difficulties involved in 
NaOH and because pH cannot 
under only 


made. 


Pecause 
accurately weighing 
be directly n easured 
1 few determinations were 

Lffeet of caleium pre sent as an umpurily During 
preliminary experiments on this method of 
synthesis, a new crystalline product was observed, 
and the same material was found in minute quan- 
tities in a considerable number of later preparations 
It consisted of slender, flexible fibers, occasionally 
2 mm or more in length. The refractive indices, 
1.584 to 1.594, were somewhat higher than those 
usually given for chrysotile, although the values 
found in the literature vary over a wide range 
Pecause of the possibility that these crystals might 
be chrysotile, a great deal of work was done in an 
attempt to increase the yield. At the same time, 
it was suspected that the fibers might prove to be 
the hydrated calcium silicate, xonotlite, because the 
refractive indices were about right for this com- 
pound. Addition of small quantities of lime to the 
magnesia appeared to have no effect. However, an 
increased vield of the fibrous crystals was obtained 
by substituting dolomitic lime for the magnesia, but 
the amount was still too small for identification 
After much futile experimentation, a product was 
obtained containing a few tenths of a milligram of 
the fibrous crystals, which could be separated me- 
chanically from the rest of the material. X-ray 
analysis proved them to be xonotlite. Although it 
was known that the magnesia contained some cal- 
cium (equivalent to 0.17 percent of CaO), the amount 
of the latter was so slight that the formation of cal- 
cium silicate crystals of conspicuous size was some- 
what surprising. It was estimated that the caleium 
must have been quantitatively converted in order to 
account for the amount of xonotlite produced. 

Miscellaneous er pe riments. The experiments dis- 
cussed above were designed to bring out in a syste- 
matic manner the effects of varving certain experi- 
mental conditions. A variety of other experiments 
were performed, in a qualitative way, in an attempt 
to promote the growth of crystals of chrysotile of 
larger size. For example, the magnesia was made 
into a paste and coated on the crucible or on other 
supports of various shapes. Again, the period of 


these conditions, 


some 


hydrothermal treatment was increased t 
Traces of various acids and of ¢ aeeaneadie of * 
metals were added, both to the silica and 
magnesia. In no case were any microsco 
visible ery stals of chrysotile produced. 

Silica in the vapor phase was found to reac: y 
magnesium carbonate at 380° C, with formato, 
chrysotile and tale. Natural serpentine was pa 
ally altered to tale by exposure to siliceous 
for 7 days at 380°, but not at 350 Under 
conditions, at 380°, olivine showed evider 
surface attack 

In another series of four experiments, an att 
was made to increase the size of fragments of n 
chrysotile asbestos fibers. Pundles of fibers we 
to a lencth of about 0.03 inch and treated 
thermally under four conditions: (a) on the surfa 
of a bed of MgO: (b) imbedded in MeO: (ce) imbedd 
in a mixture of MgO and SiO, in 3:2 ratio; (d) pack 
into cavities in magnesite. Treatment was for 
days at 420° C, with 3 ¢ of water in the bon 
other conditions were as described previously for 
vapor-transport experiments. No change in the s 
or appearance of the chrysotile fibers was obsery: 


3. Summary 


Magnesia and silica in 3:2 ratio were found 
combine under hydrothermal conditions at temper 
tures as low as 200° C. X-ray and electron mix 
scope studies showed that the product was chrysot 
but the fibrous crystals were too small to be dist 
guished as such by light microscopy. 

Chrysotile, together with tale, was also formed | 
the action of siliceous vapor on magnesia. Con 
lions governing the rate of transport of the sili 
were worked out in considerable detail. The ra 
was found to increase with temperature and _ pres 
sure and also with the initial pH. No improveme 
in erystal size resulted from this method of synthesis 

Silica was found to displace CO, from MgC‘ 
under hydrothermal conditions at 380°. 


The X-ray diffraction patterns were prepared | 
PB. M. Sullivan and G. M. Ugrinic, of the Constitu- 
tion and Microstructure Section. Their valual 
assistance is acknowledged with thanks. 
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Diclectric Properties of Teflon from Room Temperature to 
314° C and from Frequencies of 10’ to 10° c/s’ 


Paul Ehrlich? 


The dielectric constant and dissipation factor of Teflon were measured at frequencies 
from 10? to 10° c/s and from room temperature to temperatures just below the 
transition point occurring at 327° C Measurements of the d-c conductivity were made over 
the same temperature interval. Results show that the dielectric constant decreases some- 
what with increasing temperature. The relation between thermal-expansion coefficient and 
temperature coefficient of the dielectric constant is roughly as predicted by the Clausius- 
Mossotti equatior Values of dissipation factor and d-c conductivity at all temperatures 


first-order 


are very low (less than 2 10~* and 2X 10-" mho/em, respectively), and the dielectric con 

ibedd stant is independent of frequency at all temperatures. 
| pack 
s [Or 
» bon 1. Introduction 3. Methods 

for t 

the s lielectric behavior of Teflon, an organic high 3.1. Electrical Measurements 
bserve with outstanding temperature stability, is 


zated over a wide temperature range with 
ent especially designed to permit measure- 
at elevated temperatures. Previous qualita- 
vestigations |1]*° have indicated the absence of 


The measurements of dielectric constant and dis- 
sipation factor were carried out by means of a sub- 
stitution method, changes in sample capacitance with 
temperature being followed by balancing with the 


we A ant dielectric losses in Teflon up to its first- internal bridge capacitor connected in parallel with 
. micro. transition point at 327° C No quantitative the sample holder. 
rvsot tigation has been made of the variations of die- In the experiments to be described, the sample 
» diati: constant and dielectric loss with temperature Capacitance at room temperature was found by mak- 
requires the use of a loss-free micrometer ‘Ng Measurements with the sample between the elec- 
med | le holder having the characteristics discussed trodes and with the sample removed from the elec- 
Cond where [2]. trodes. After the temperature was raised, the 
16 silir lhe following symbols are used in this paper. sample was kept between the electrodes at all times, 
‘he rat subscript 0 refers to room-temperature values and drifts in the bridge were compensated for by 
id. pres disconnecting the sample holder after each reading 
ines and by reestablishing the balance with the internal 


nthesis 


MeC0) 


ared | 


onstitu- 


raluabl 


2. Symbols 


the dielectric constant. 

the dissipation factor. 

tan 6=the loss factor. 

frequency in cycles per second. 

2rf. 

relaxation time 

the capacitance of; a sample of air of the same 
dimensions as the dielectric measured. 

the electrode diameter (2 inches). 

the sample diameter 

the sample thickness 

the capacitance reading on the internal bridgs 
capacitor with the sample holder connected. 

the capacitance reading on the internal bridge 

capacitor with the sample holder disconnected 

Cg—C computed at any constant temperature. 

the difference between the room temperature 
value of Cp and that obtainea at another tem- 
perature. 

the computed value of the capacitance of the ring 
of air located between sample and electrode 
edges, neglecting edge effects. 

the difference between the value of C, at any 
temperature and the room temperature value. 


| by Office of the Chief of Ordnance, Department of the Army. 
idress, Department of Chemistry, Harvard University, Cambridge, 


n brackets indicate the literature references at_the%end of this}paper. 


bridge capacitor. The disconnecting was accom- 
plished by means of the special switch built into the 
sample holder. Details of the measuring equipment 
are described elsewhere [2] 

If it is assumed (1) that the sample expands uni- 
formly in all directions, and (2) that the edge effect 
at all temperatures is effectively eliminated, we 
obtain 


O.OSS544r 


AC yw 4s [Pra 


and the dielectric constant at temperature is 


found to be 


any 


8,/8 [eg +(ACn+AC,)/C 


J. 


3.2. Sensitivity and Accuracy of Measurements 


The sensitivity of the dielectric constant and dissi- 
pation-factor measurements was determined by the 
sensitivity and stability of the bridge elements, as 
well as by the ratio of sample capacitance to that of 
the entire arm, consisting of leads, sample holder, 
sample, and substitution capacitor. It is estimated 
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that, for the experiments cited, the overall accuracy 
of the dielectrie-constant computations is probably 
not determined by the accuracy and sensitivity of 
the dimensional measurements made at room tem- 
perature or by the bridge balance, which ceuld be read 
to +0.002uuf. The overall accuracy is determined 
bv deviations from the assumptions made, 
namely, that warping does not take place, that the 
sample expands uniformly, and that the edge correc- 
tions at all temperatures are effectively eliminated, 
Small deviations from the first two of these assump- 
tions, which are not readily recognizable, are expected 
to exceed the sensitivity of a room-temperature 
measurement. It is therefore difficult to make a 
reliable estimate of the accuracy of the dielectric- 
constant measurements carried out at higher temper- 
atures. As shown later, the agreement between the 
several runs was very satisfactory, but an answer to 
the question of accuracy has to wait for further exper- 
iments on samples of widely varying dimensions. 
The best indication that the dielectric-constant 
values computed cannot be far from correct is the 
agreement between the observed and predicted re- 
lationship between thermal-expansion coefficient and 
the temperature coefficient of the dielectric constant. 
The estimated accuracy of the dissipation-factor 
1x10 The d-c conductivity 

smallest observable value of 


basic 


measurements 1s 
never exceeded the 
210 mho/em. 


3.3. Preparation of Specimens 


Molded Teflon undergoes drastic dimensional 
changes when reheated. It was expected that dimen- 
sional stability could be attained by conditioning the 
samples at a temperature above its 327° C first-order 
transition point. Complete stress relaxation and 
disappearance of all crystallinity at this temperature 
has been reported [1], and it was therefore expected 
that the conditioning process would relieve internal 
strains introduced during molding. This was borne 
out by the fact that samples that underwent very 
drastic changes during conditioning for several hours 
at 350° C exhibited no further obvious changes when 
the process was repeated and did not appear to warp 
during the electrical measurements. Prior to ma- 
chining, the samples were therefore kept at 350° C 
for several hours. They were then machined in a 
vacuum chuck to give disks approximately 1‘ in. in 
diameter and approximately \ in. thick. As meas- 
ured by calipers, the faces were flat and parallel to 
+0.0003 em and the diameters constant to +0.005 
em, the latter value being determined by the sensi- 
tivity of the caliper. 

Previous experiments had shown that for accurate 
measurements over a temperature range it was 
essential to apply secondary electrodes to the samples. 
Of the materials tested, evaporated gold was found 
to be the most satisfactory electrode material, and 
this was used in all the experiments to be cited. 


4. Results and Discussion 


Measurements were made on samples A anc B 
from different sections of a molded rod. Sample | 
was reannealed as previously described to give sam 
C. Sample D had been machined from 
molded separately and came from a different lot ¢} 
the rod from which the other samples were cut. T 
values of the dielectric constant at any tempera 
were not frequency dependent, as will be demoy. 
strated presently, and the values measured at ro 
temperature for samples A, B, C, and D were 2.0 
1.994, 2.016, and 2.022, respectively. The densit 
the samples was 2.13. 

Table 1 gives results for a typical run (samp) 
carried out over the entire temperature and freque; 
ranges. Values of the dielectric constant have by 
calculated with the aid of eq (1) and (2). It is 
that the values of the dissipation factor, which , 
of the order of the sensitivity of the bridge, do ) 
show a trend, either with frequency or with temp 
ture. The dissipation factor can, therefore, be rep 
sented as having the value of (14£1)*10~* over ¢/ 
temperature and frequency ranges investigat 
Inspection of table 1 also indicates that there is | 
significant variation of the dielectric constant 
frequency at any one temperature. Such a variati 
would not be expected in view of the very low valus 
of the dissipation factor. The frequency depender 


of the dielectric constant can be predicted from t| 
observed values of the loss factor by extending 
procedure already described [4] to materials wy 
frequency-independent losses, giving 


de’ 
=2.303Ae", 
d log f 303 Ac 


where 


whence 


de’ 
== 1.476". 
d log f bs 


Thus, for Teflon (e’=2, tan 6=1X10~*), a chang 
in dielectric constant of 0.001 over the frequen 
interval from 10? to 10° c/s is computed. This cor 
firms the conclusion, already surmised from tabl 
that there is no significant variation of dielectr 
constant with frequency at any one temperature 
One would expect, therefore, that the variatior 
dielectric constant with temperature would be caus 
essentially by the variation in the bulk density of | 
material. It becomes clear that the accuracy of | 
dielectric-constant computations depends critical) 
on the validity of the assumptions underlying eq 
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an accurate determination of the sample 
ess at all temperatures. Figure 1 shows a plot 
ve sample thickness versus temperature. It 


en that the reproducibility is satisfactory and 
the experimental values agree with those reported 


of the manufacturers of Teflon. This may be 
as evidence that warping of the sample did 


take place to an extent that would seriously in- 
nee the results. 


ire 2 is%a plot of dielectric constant versus 
iture. An average value of the dielectric 
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constant is used, valid over the entire frequency range- 
and room temperature values of 2.000 have been as 
sumed for each sample in order to obtain a more 
accurate idea of the reproducibility of the temperature 
runs. The values used in the plot have therefore been 
obtained by merely adding or subtracting a small con- 
stant correction from the experimental points at all 
temperatures. As in figure 1, the reproducibility is 
satisfactory 

The Clausius-Mossotti equation predicts a relation 
between the temperature coefficient of the static 
dielectric constant and the linear thermal-expansion 
coefficient, a, which, in the case of a nonpolar ma- 
terial like Teflon, one might expect to be applicable 
to the dielectric constant at audio and radio fre- 
quencies. Similar relations can_be derived for other 
theories, but we shall assume the validity of the 
Lorentz field and therefore of the Clausius-Mossotti 
equation, obtaining [5] 


de’ 
dT 


In the case of Teflon, the value of de’ /AT)/ a; pre- 
dicted on the basis of the above equation would 
therefore range from roughly 3.85 at 100° C. to 
roughly 3.5 at 250° C. The experimental values 
range from 3.6 at 100° C to 3.0 at 250° C 

In view of the assumptions underlying the evalua- 
tion of the experimental dielectric constants, this 
comparison would seem at least not to be in disagree- 
ment with the Clausius-Mossotti theory, although 
a truly quantitative test on this basis would seem to 
be impossible. 


The sensitivity of the bridge is much too low to 
obtain even a qualitative idea of the behavior of the 
dielectric losses as a function of temperature and 


frequency. Table 1 indicates that the values of the 
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Figure 2 Dielectric constant of Teflon at any frequency 
hetween 102 and 10° c/s as a function of te m perature 
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dissipation factor are of the order of 10 No satis- 


factory data on nonpolar polymers with losses of 
that order of magnitude have as vet been published 
and theories [6, 7] put forward to explain the nature 
of such losses have therefore remained untested 
Quantitative information on phase transitions im 
Teflon ts available both on the basis of dilatometric 
8] and of thermodynamic measurements [9] 


These 
investigations agree as to the occurrence in the 
room-temperature range of a double first-order tran- 
sition, both steps of which presumably occur between 
crystalline forms, the less pronounced step occurring 
at about 30° C, whereas the one involving the larger 
volume change occurs at about 20° C The volume 
changes associated with this transition were not imdi- 
cated by the thickness measurements made on the 
sample in the holder because these measurements were 
very crude compared to dilatometric measurements. 
The fact that there appears to be no drastic change in 
dielectric constant in this appear 
surprising if one assumes that the transition occurs 
The high- 


can be 


region does not 


between essentially nonpolar crystals 
temperature melting 
easily indicated by the dielectriv 
is associated with dimensional changes of sufficient 
magnitude to result in large capacitance changes 
There was considerable warping of the sample at 
this temperature, and attempts to carry out repro- 
ducible measurements through this transition range 


transition at 62, & 
measurements as it 


were abandoned. 


5. Summary 


The variation of the dielectric properties of Tos 
from 24° to 314° C can be accounted for in terms 
the variation of the bulk density of a nm 
material, 


The author expresses his appreciation to A. jj 
Scott for helpful advice throughout this invest 
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E of Solutions of Potassium Tetroxalate from O° to 60° C 


Vincent E. Bower, Roger G. Bates, and Edgar R. Smith 


The National Bureau of Standards conventional activity scale of pH is defined in terms 
of a series of standard buffer solutions prepared from certified materials issued as NBS 
Standard Samples In order to increase the accuracy of measurements at low and hig! 
pH, standards of th acidity and high alkalinity are needed to supplement the four now 
available, which cover adequately the intermediate range from 3.5 to 9.2 

Potassium tetroxalate dihydrate, chosen as the fifth pH standard, is a stable, crystalline 
substance composed, in addition to water, of equal molar amounts of oxalic acid and potas- 
sium acid oxalate Its solutions are therefore well buffered in the concentration range 0.01 
to 0.1 M and are highly acidic, with pH ranging from 1.5 to 2.2 Eleectromotive-force data 
were obtained for 12 mixtures of potassium tetroxalate and potassium chloride at 13 tem 
peratures from 0° to 60° C, and pH values were assigned to 0.1-, 0.05-, 0.025-, and J0.01--M 
solutions of potassium tetroxalate without added chlorid The 0.05-M solution, which has 
a pH value of 1.68 from 15° to 35° C, is recommended as a standard for meas 


low pH 


l. Introduction In eq (2), Fand Es represent the electromotive force 
of cell 1 when the unknown and the standard solution 
are two primary reasons for the establish- whose pH is pHs), respectively, are in contact with 
pH standards. The first and more impor- _ the electrodes 
to provide a reference for the adjustment of The hydrogen and hydroxy! ions have much greater 
equipment to a common basis, making mobilities than other ions. Furthermore, the mobil- 
true reproducibility of pH measurements ities of these other ions are all of the same ordet 
me to time and from place to place. The For this reason, the liquid-junction potentials depart 
purpose of pH standardization is to simplify appreciably below pH 3 and above pH 11 from the 
erpretation of the measured pH and to make — relatively constant values maintained in the region 
niform over the entire range of the scale. of intermediate pH At high acidities and high 
first of these objectives can be met by the  alkalinities, therefore, the measured pH may not be 
sal adoption of a single pH standard. The susceptible to the same simple interpretation given 
s not so easily achieved. The unknown and it near pH 7. When measurements below pH 3 and 
ile potential at the liquid junction, indicated above pH 11 are to be made, this defect can be 
vertical line, in the practical pH cell remedied to some degree by the choice of standards 
that are themselves highly acidic or highly alkaline 
H, (or glass electrode), solution X! 3.5 M or (1) 
saturated KCI. He.Cl.: He | To permit as far as possible a consistent interpre- 
a tation of the pH over the entire scale, the National 
the interpretation of measured pH values in Bureau of Standards is establishing a series of stand- 
f concentrations or activities very difficult, ards extending from low to high pH. ~The standard 
it is justifiable to assume that the practical pH is based upon emf measurements of cells free from 
son the seale defined by the standards onlv @ liquid junction and its attendant uncertainties 
, xf The standard pH scale is thus defined in terms of 


solution X is a dilute aqueous solution of 
several fixed points in much the same manner as the 


e jonized solutes. : om 

\s is well known, the liquid-junction potential re- International remperature Seale ; : 
is from transfer of the positive and negative ions The pH of the standards is defined formally by 
ransference numbers of all the ions are differ- pHs log au, 3) 

ind possibly also the concentration gradients of 
several ions through the boundary. the net charge Where @q is a conventional hydrogen-ion activity. 
nsferred is rarely zero It is impossible to make Also by definition, @y fyi, n which fy and my are 
correction for the variable liquid-junction the activity coefficient and molality of hydrogen ion 
ll. Hence, an unchanged liquid, junction Inasmuch as the activity coefficient of a single ionic 
il on substitution of the standard (S) for the Species cannot be determined, the numerical values 
n solution in the cell is implied in the opera- of pHs must have a conventional or defined basis 
inition of the practical pH value: The NBS pH values rest upon the average of several 
reasonable estimates of the individual activity co- 
E—E.)F efficient. At ionic strengths less than 0.1, the pHs 
pH=pHs+ >» . 2) eid 

RT |\n 10 Figures in brackets indicate the literature 1 ene tt end t} reat 
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values computed in these several ways are the same 
within about +0.01 unit [2, 3, 4] 

Potassium tetroxalate dihydrate, a stable acid salt, 
is a promising standard for the acid region. Its 
aqueous solutions, formed by weighing a single solid 
substance, contain equal molar amounts of oxalic 
acid and potassium binoxalate and are therefore well 
buffered. The pH is of the same order of magnitude 
as —log A), where A, is the first dissociation constant 
of oxalic acid This acid group is moderately strong, 
with A,=0.0536 at 25°, according to Darken [5]; 
hence, solutions of potassium tetroxalate with pH 
less than 2 can be prepared 


2. Method 


The method by which the standard pH values 
pHs) were assigned to solutions of potassium tetroxa- 
late has earlier” There are three 
steps in the procedure. First, the emf, /, of cells of 
the type 


been described 


tetroxalate (” KCI (m, 


AgCl; Ag 


atm), Kx 


Pt; H, (g, 1 


was measured and, with a knowledge of m., the 
molality of chloride ion in the buffer solution, values 
of the quantity pwH calculated by the equation 

fa foumn . o°) F/(2.3026RT 


log Mo, 


pwH = —log 


where £° is the standard potential of cell 4 and the 
other symbols have their usual significance 

The second step was to determine pwH?®, the limit 
approached by pwH for each buffer solution as the 
molality of added chloride was reduced. This was 
accomplished by the determination of pwH for 
chloride molalities of 0.015, 0.01, and 0.005, plotting 
pwH as a function of Mey that is, Me and extra- 
polating to mo,=0 

The third step was to compute pH, from pwH 
From the definitions of pHs and pwH, we have! 

pHs=pwH°+loz fa, (6) 

where fc is the limit approached by the activity 
coefficient of chloride ion as the composition of the 
buffer-chloride mixture appgoaches that of the pure 
buffer solution. A nonthermodynamic assumption 
must be introduced at this point to permit the ionic 
activity coefficient to be evaluated. 

The several reasonable conventional definitions of 
this coefficient that have been employed earlier for 
the computation of pHs [6, 10, 11, 12] are all con- 


? See, for exemple, the assignment of pHs values to solutions of potassium hydro- 
gen tartrate [6] 

' The values of F 
of the standard potential were der 
Ehlers [8] as recalculated by Hamer 


and 2.W26F T’ Fin absolute volts are listed in [7 The values 
ved from the measurements of Harned and 
and Acree [9 


Burton, 


sistent with the” Debye-Hiickel equation 


—log fo Avu 
1+ Ba vu 


in which w is the ionic strength, a; is the 
parameter, and A and B are constants for a 
solutions at a particular temperature. Th 
sumptions differ essentially only in the valu 
and the manner in which this value was 
By combination of eqs (6) and (7 

pHs pwH —Ayu l t Ba ve 
From an examination of the mean activity 
cients of hydrochloric acid and chloride salts 
be concluded that a, for chloride ion probab 
between 3 and 6. When the concentration 
tetroxalate is 0.05 .VW/, this uncertainty (+1.5 
amounts to less than 0.01 unit in the assigned pl 


3. Experimental Procedures and Results 


The potassium tetroxalate dihyd: 
KH,(C,O,).-2H,O) was a composite sample of 
250-lb lot purchased for issuance as NBS Stand 
Sample 189. The material was slightly m 
Before it was dried it was found 99. 
99.7 percent by titration with a 
0.17— solution of sodium hydroxide to th 
point of phenolphthalein. After drving three s 
ples for 2 to 3 hours in air at 50° C, they assay 
99.93 + 0.02 percent. The salt should not be hea 
above 60°C [1]. The results of other 
tests [14] were as follows: 


to assay 
carbonate-{ 


stand 


Chloride: Less than 0.002 percent 


Sodium: Less than 0.02 percent 

Calcium: Less than 0.005 percent 

Iron: Less than 0.001 percent. 

Ammonia: Less than 0.002 percent 

Sulfur compounds: None perceptible 

Heavy metals: Less than 0.002 percent (as P! 
Substances darkening sulfuric acid: None 

The potassium chloride was a bromide-free fus 
sample prepared in the manner described earlier 
The specific conductance of the water used was 
the average, 0.810-° ohm='em™. Dissolved 
was removed from the solutions by passing nitrog 
through them for 2 hours. 

Four concentrations of potassium  tetroxala! 
namely 0.1, 0.05, 0.025, and 0.01 m/, were study 
The emf values, corrected in the usual way & 
partial pressure of 1 atm of hydrogen, are summari 
in table 1. Each recorded figure represents 
average for two pairs of electrodes in the same ‘ 
These potentials are in good agreement with sim 
but less extensive measurements made in an eat 
investigation [1]. 

‘ Values of the Debye-Hiickel const»nts from 0° to 100° C will be f 
In the formula for * (bottom of p. 1765), “10-*" should read “10 





TABLE 1. Electromotive -force values for cell 4 from G° to 60° C. 


14004 
44009 


>sults 
Val 
le Oo; {nnd 
stand 37 12 54 3 5 4048 
f anal i268 it 
it - ~~ 
$6952 723 4778 
OU} | OT re wt 4h900 712 7237 T3225 7 r ; wl HS 4773 
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late-! ‘7 
; ‘771 
he er 323 7 
ree 8 j ‘ 4 48670 3 O75 112 2 5 W557 wei 
" Ose . 4 , s4Hu 


. iso 
assay 
e heat 


stand 
values of pwH°® obtained are summarized in table 2 


The accuracy of pwH?® is judged to be about + 0.003 
5° C and +0.005 at 60° C 


4. Determination of pwH° 


pwH was calculated by eq (5) from each emf at 0° and 
alue recorded in table 1. In order to obtain pwH 
the value of pwH for each tetroxalate solution in the OLALITY OF Ci/u 
absence of chloride, two types of extrapolation, ’ 0.10 
llustrated in figure 1, were tried. The pwH values . T 
the 0.1-m solution of potassium tetroxalate were 
first plotted as a function of the molality of chloride 
\ straight line, the lower line of figure 1, was drawn 
ind extended to zero molality of chloride. The 
line was drawn through the pwH values for 
lutions containing the three different amounts of 
chloride but of lower tetroxalate content, so that the 
-oxalat three points correspond to a constant ionic strength 
stud equal to that of the chloride-free buffer solution. 
av te Th pwH, obtained by interpolation, is plotted a 
mari function of me/u, the fractional contribution of the 
nts potassium chloride to the ionic strength. 
me ci The two types of extrapolation gave the same 
result, within the error of the determination. Hence, 
the simpler and more convenient treatment, illus- 
trated by the lower curve, was used. The plots of 
pwH asa function of the molality of added potassium 
chloride at 0°, 25°, and 60° C are shown in figures 
-, 3, and 4, respectively. The positions of the lines pyoure 1. Extrapolation of pwH tofzero concentration 
were determined by the method of least squares. The potassium chloride by two different methods 





was 
lved 
nitroc 
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Fiaure 2 
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FIcuR! pwH of solutions of potassium tetrora 


inctre f the molality of added potas im chioride at ¢ 








5. Calculation of pH, 


The ionic strength of a solution of potass 
tetroxalate of molality m, is given by 


vi ny Ng. 








The dissociation of the second acid group of oxa 
acid in tetroxalate solutions is so slight that 

contribution need not be considered [16]. Tl 
hydrogen-ion concentration, my, was obtained w 
sufficient accuracy by the method of success 
approximations directly from pwH° with activit 
coefficients estimated by the Debve-Hiickel equati 

















log my=pwH°?+log (fufer pwH 


pwlt of 
i the mo 


H° f , ' The values of the ionic strength at 0°, 25°, and 60° | 
nw | d 0 of? 128 ; : 4 
;, obtained in this manner are given in table 3 


TaRLE 3. Jonic stre naths of four solutior x of potas tu 
tetroralate 





buffer substance, although the change of pH upon 
dilution increases near the ends of the pH scale 
It is possible, therefore, to take advantage of the 
convenience of volumetric methods for the prepara- 
tion of buffer solutions. It is evident from figure 5 
that the pH of tetroxalate solutions at 25° changes 
about 0.2 unit when the molality drops from 0.05 to 
0.025. Two determinations with a 50-ml picnometer 
showed the density of the 0.05-molal solution to be 
1.00320 at 25° C: hence, the 0.05-molal (m) solution 
is 0.0496 molar (17). This difference of 0.8 percent 
in the tetroxalate concentration between 0.05-m 
and 0.05-\/ solutions corresponds to O0.OOL to 0.002 
pH unit, an entirely negligible amount in comparison 
with the uncertainty in the assignment of pHg 
It is therefore justifiable to assign the pH, values 
determined as described above for 0.01-, 0.025-, 
0.05-, and 0.1-m solutions also to the solutions of 
numerically the same molarity, as has been done in 


table 4. 








1 1 l 
0.02 0.04 0.06 0.08 
MOLALITY OF K TETROXALATE 





6. Tetroxalate Buffer Standards 


The 0.05-M solution is recommended for use as a 
standard for pH measurements in the acid range 
In common with other strongly acidic solutions, it 
has a high buffer value and a high dilution value. 
In other words, its pH is little influenced by chance 

Hs at 25°, caleulated by eq (8) for three contamination with acid or base but is relatively 
ola namely, 2, 4, and 6, is plotted in figure sensitive, in comparison with buffer solutions of 


pH of solution 7 potassium tetroralate 
as a tune 


tpH 


curves converge as the ionic strength moderate acidity, to dilution or concentration of the 
ses. Hence, it is possible, for molalities of buffer components. Nevertheless, the dilution value 
ilate not exceeding 0.05, to select ‘an average 18 not high enough to impair the usefulness of tetroxa- 
it will represent the values on all three curves late buffer solutions in any significant degree. 

O.Ol unit Nevertheless, when the mean The Van Slvke buffer value, 8, is defined as the 
coeflicients of hvdrochloric acid in its pure differential ratio dB dpH [18] and represents the 
s solutions and in solutions with added potas- number of equivalents of strong base (/3) required, 
Its are fitted to equations of type (7), the without volume change, to increase the pH of 1 liter 
parameter is usually found to lie between of a particular solution by | pH unit. Approximate 
~ For this reason, the pHs given in table 4 buffer values are obtained experimentally by measur- 
verage of those computed by eq (8) with ing the pH change produced by the addition of small 
ind a;=6. In view of the arbitrary selection Volumes of solution, each of which contains APB 
in uncertainty of +0.01 is assigned to the equivalents of concentrated, strong alkali In this 
for 0.01-, 0.025-, and 0.05-M solutions and way, the buffer value of 0.05-\/ potassium tetroxa- 
to those for the 0.1-/ solution. late at 25° C was found to be 
pH of buffer solutions is not normally very 

to changes in the total concentration of dB AB 0.07 11) 

dp ApH rte 
pHs of 
The pH change was determined with the glass elec- 
trode. 

The dilution effect is conveniently expressed as 
the increase of pH on dilution of the buffer solution 
with an equal volume of pure water. For the 
0.05-M solution of potassium tetroxalate, this 
quantity, designated ApH,» is 


ApH, 2 0.19. 

















pH of 


Fieure 6 


The buffer value of the 0.05-\/ solution of potassium 
acid phthalate is about 0.02, and its dilution value is 
+0.05; hence, the pH of the tetroxalate solution is 
altered less than that of the phthalate solution by 
addition of acid or alkali but is more sensitive to 
dilution. 

As shown in figure 6, the temperature coefficient 
of pH is small and positive for all four tetroxalate 
solutions in the range of the measurements. The 
pH of the 0.05-M solution is 1.68+0.01 at all tem- 
peratures between 15° and 35° C 

A solution of potassium tetroxalate saturated at 
at 0° C was used by Schomaker and Brown [19] in 
their standard quinhydrone half-cell. Hitchcock 
and Taylor [20] selected the 0.1--/ solution of this 
salt as a useful pH standard and assigned to it pH 
values of 1.480 at 25° and 1.495 at 38° C. These 
pH values are 0.040 unit lower than those given 
in table 4, or about two and one-half times the 
estimated uncertainty of the latter. This difference 
corresponds to 2.4 my. 

The determination of Hitchcock and Taylor was 
based on the emf of cells with a liquid junction, and 
the primary reason for the difference between the 
values of the two investigations is probably the 
residual liquid-junction potential, which becomes 
appreciable at low pH values. In other words, 
the potential across the boundary between the 


bridge solution (saturated potassium chlori 
0.1-\f potassium tetroxalate differs fro: 
between the bridge and other buffer solutions why, 
pH values lie closer to neutrality. 

Potassium tetroxalate dihydrate can be drie 
air at l-atm pressure at temperatures up to aly 
60° C without loss of crystal water. The 6.95.1 
solution is prepared by dissolving 12.70 g (air weig| 
in water and diluting to 1 liter. In an experim 
designed to reveal the effect of drying the salt 
solutions were prepared by adding the same wei 
of wet and dry salt to equal amounts of a soly 
of potassium chloride. Although the wet san 
contained about 0.5 percent of water, the emf 
cells prepared from the two solutions displayed 
significant differences. Hence, it appears that dryig 
is usually unnecessary. Textroxalate solutions 
satisfactorily stable without the addition of prese: 
tives. In order to minimize the effects of evap 
tion and accidental contamination, however 
recommended that the standard be replaced ey 
2 months. 
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On Small Disturbances of Plane Couette Flow 


Wolfgang Wasow’ 


The Orr-Sommerfeld equation for the function ¢(y), which appears in the component 
¢(y)exp|ia(x—ct)| of the stream function, is analyzed in detail for the case of plane Couette 
flow. A set of solutions is found whose asymptotic behavior as aR—o can be calculated 
in the whole complex y-plane (FR is the Reynolds number A disturbance with given a is 

’ Welg proved stable, for any a, if aR is sufficiently large In addition, an investigation is made 
of the asymptotic properties of an allied equation for which, in contradistinction to the 
singular Couette case, an inner friction laver does exist 


solut 
t san 
e emi 
laved l. Introduction stable, if one has Im ¢<0 for all possible eigenvalues. 
at drving Otherwise the flow is unstable. 
tions This approach has been very successful for a 
poeae me questions connected with the stability of number of velocity profiles. It has led to critical 
evan ; evnoidads numbers mt are in good agreement with 
vay athe Gon ' Reynold I that 1 ay t witl 
rer, Us ‘e betw two infinit ‘allel pl experimental evidence, [4, 5, 6]. A mathematically 
ed ev ' Sa 0865S “ . = fill yh et Planes difficult part of this method is the integration of the 
aed fund - lif. - 0. wit ; . a hile differential equation (2). For large values of aR 
esslDle ule , and lu one plane is at rest Whli¢ asymptotic methods can be used whose validity 
her moves with constant velocity, then the has been established in [7, 8 and 10} . 
vent yo’ ——— with — Navier ; a Couette flow, however, although it has the simplest 
1ations i : caites ae cies eclaaier a file. velocity profile, plays an exceptional role in these 
nary parallel Moto 1 lear Vvelocit: rome. PK ie ‘re “§ . s of [7 ‘ 
plane Couette flow. Experiments indicate theories. First, the proof of [7, § and 10} do not 
oy that this fi “is , : alle ie ’ cover this case. Second, all investigations seem to 
= ak = ig ay “ge le tor Thais indicate that Couette flow is stable, in the sense 
org Ss ? , Fy Ss ‘ , » oO “o 
“sega ae pte, Tat 1s, Turbulent, lor Mgh specified above, at all Reynolds numbers, although 
= no mathematically complete proof of this has been 
oldest and best developed theory of hydro- given (ef. [2], p. 256). Third, in contrast to other 
stability for plane parallel flows, not only profiles, no so-called “inner friction layers’’ 3, 8] 
Couette flow, studies the stability of disturbances oeeur in the theorv of Couette flow. 
» sma that they can be regarded as solutions of a This paper will consist of three parts, correspond- 
ized differential equation. Itisfurtherassumed jing to these three points. 
stream functions of such disturbances can (a) The first part (section 2) is in the nature of a 


irpose of this paper is to analyze mathemat- 


“ined by superposition of functions of the form supplement to [9] and [10] A fundamental svstem 
of the differential equation (2) for the case of Couette 
d(yjee?-™, (1) flow is constructed whose asymptotic behavior, 
as aR->o, can be calculated in the whole complex 
the main flow is taken to be parallel to the z-axis. y-plane. 
y is positive, ¢c and ¢(y) are complex, and ¢ is b) In the second part (section 3) the eigenvalues 
The references [1, 2] contain fuller ac- e¢ of (2), corresponding to possible disturbances of 
ts of this theory. It reduces the stability prob- — the flow, are studied as functions of a and aR. It is 
; eigenvalue problems for the differential proved that yor any given a a disturbance of type (1) 
h NBs %, juation of Couette flow is stable, if aR is sufficiently large. 
h NBS If the physically plausible conjecture is accepted, 
a'¢—iaR| (w—c) (¢’’ — a’) —w’’¢]=0, that above a certain frequency a there is always 
(9) stability (cf. [3], part ILI, p. 284), the stronger 
‘i statement can be made that, uniformly in a, there 
is a lower bound for aR above which the flow is 
stable. 

On the other hand there are the inequalities of 
Synge ([2], p. 258), which imply the stability of 
Couette flow for sufficiently small values of R. 
Even if taken together, these results do not, of course, 
constitute a complete mathematical proof of stability. 
But the author believes that they go farther towards 
this goal than previous efforts. 

C. Morawetz has shown [13] that the possible 

eigenvalues ¢ of the stability problem fall asymp- 
ireau of Standards, Los Angeles, Calif . a , Ns , 
ekets indicate the literature references at the end of this paper. totically into three classes: (1) Eigenvalues that 
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er w(y) is the velocity of the main flow and R 
enotes a Reynolds number of the flow. 
A solution (eigenfunction) of (2) that satisfies the 
omogeneous boundary conditions of a particular 
ow problem under consideration yields a stable or 
disturbance of that flow according as the 
nding eigenvalue c has a negative or positive 
part. For a given value of R the flow is 


tion of this paper was sponsored (in part) by the Office of Naval 





approach those of the limiting inviscid problem as 
ako. The corresponding eigenfunctions approxt- 
mate those of the inviscid problem except possibly 
2) An iniinity of 
close to a certain 

They are in no 
equation 


in certain “friction layers.” 
eigenvalues lying increasingly 
curve in the e-plane as af 

wav related to the inviscid 

The corresponding disturbances are always stable 
Eigenvalues that approach the points w(y,) or 
wy), as aR—-o, where y, and y are the ordinates 
of the flow boundaries. The actual existence or 
stabilitv of eigenvalues of this third kind was not 


differential 


discussed in [133] 

Since no eigenvalues of the kind 
plane Couette flow, its asymptotic stability, in the 
sense stated above, is proved if the disturbances 
of the third kind are shown to be stable for suffi- 
ciently large values of R. This will be done in 
section 3 

c) In all cases examined so far, the actually pos- 
sible disturbances of laminar flows are approxima- 
tions to solutions of the inviscid differential equation, 
except, in general, in certain sectors of the complex 
y-plane, where they diverge exponentially, as aR— « 
ef. [8, 9, 10] The results of section 2 show that 
no such sectors of divergence exist in the Couette 


first exist’ for 


cause 

In section 4, the mathematical reason for this 
absence is analyzed. Although not immediately 
applicable to physical problems, the results of that 
section may turn out to be of more than purely 
mathematical interest. 


2. Solution of the Differential Equation 


In terms of the dimensionless quantities used in 
(2), the velocity profile of Couette flow is, for an 


appropriate deiinition of R, 
uw Y 


The further substitutions 


oly u(z 


change (2) into the simpler form 


— ae 
u 2aru 


tatu +22 (u"’ au 0, (7) 
The theory of [10] does not apply to this differential 
equation, because condition (3.3) of that paper is 
not satisfied. But (7) can be treated directly by 
the method of Laplace transformation used in [9}. 

By a straightforward application of Laplace’s 
method, it is seen that contour integrals of the form 


are solutions of (7), provided C is a path that is 
either closed or that extends to infinity in such a 





Five such paths ar 
where the sectors in wi 
Re (~*f) >0 have been shaded. The paths ( 
(j= 1,2,3) are the same as in [9 and 10]. The » 
tions corresponding to these contours will be ea 
A,(2,d), U(2,) and V(z,A). Clearly, 


A,+A,+A V. 


way that Re(“f) <0. 
cated in the figure 1, 


Let S denote any bounded neighborhood of 
origin in the z-plane. Then a short calculat 
shows that, in S, 


cosh az 


‘ . 
sinh, az. 
a 


The asymptotic calculation of the solutions A 
proceeds exactly as in [9]. It is amply suffi 
therefore, to restate those results. 

Let — and » be defined by 


9 


— \(—2 


=3 


2,3) the three rays Re § 


and denote by C, (k=1 
numbered so that 
4k 2 


arg < ™—3 arg A, on C;,. 


The closed sector of S bounded by C,,(j#k 
of central angle 27/3 may be called S,. If & is 
arbitrary positive constant, one has then 


(a) For || >&, z in S—C,, 


A, (2,4) = er? (: 
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determination of — for which Re ¢& 


ESO0inS 
0 in S—S, must be taken 


Ld 


(*(A,) is obtained from C(. by a rotation 
i 2/3 arg X 
Near C, the asymptotic behavior of A, for 
follows immediately from (9), (11) and (14 
se asymptotic relations may be indefinitely 
y differentiated. Without much additional 
can be shown that (14) and (15) are first 
‘imations of infinite asymptotic series progress- 
powers of &' and \~**, respectively. 





Figure 2 


angle of €-° in formula (14) is 2/2—2xk/3 
ryat the points of the angular bisector of S— ¢ ‘be 
other points of the domain S—C,, it is then 


A convenient fundamental system for the present 
discussion is formed by the functions A;, Ay, UU, V. 
Then the eigenvalues are the solutions of the equation 


ed by analytic continuation. The angle of 


s —2/3 arg \. The angle of \ will be fixed in 
ordance with (5) by setting 


3. Stability Problem 


coordinate system and the definition of the 
nolds number can, and will, be chosen so that 


Two ranges of ¢ have to be distinguished according 
oundaries of the flow are at 


as n.=A*"e remains bounded or is unbounded as 
Ao. Ifn,is very large, the asymptotic expressions 
can be used at both end points, and arguments 
analogous to those employed in [13] show that the 
corresponding value of ¢ cannot have positive im- 
aginary part. 


Y:=VU, Ys , 17 (14 


boundary conditions for the disturbance are 


s, p 134 


If ». remains bounded as \->@ the expressions 
(15) are to be used at z=1—e and the expressions 
(14) for z c, that is, n=n Setting, for abbrevia- 
tion 


(0 o(1 o' (0 d’ (1 0. (18 


ms of the variables introduced in (4), (5) and 
this becomes 


, aetnt eBay 
ull ° um’ | 7 , 0. (19 ‘ 


The differential equation (7) with the boundary 
ondition (19) constitutes the eigenvalue problem of 

plane Couette flow. In view of the results of [13] 

lescribed in section 1] only those eigenvalues e=c(a,\) 
ed to be studied, for which either 


one finds, after a short calculation, that F' 
becomes, in this case, 


9i(ne) +X? Ea, d,€,ne) =0 (24) 


lim c(a.d) =0 where E(a,d,¢,n-) is bounded for fixed a, as X >, 
a : provided », remains bounded. Using Rouché’s 

theorem, as in [13], it can be proved that, if 7= 1 is 
lim ¢(a,d . a root of the equation 


A~ 


, 5 eee gi(n) =90, 
notation differs from that of [13] in that there ; 


j 


: then 
and that the names of the sectors S,; and 


nterchanged.) The case (20) will be con- ; 
first. Using (13) and (16), the sectors S, are that is, 
have the appearance indicated by the figure 2. c(a,d) 


+ Or 
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If c(a.d)—>1, the calculations are analogous. The Here the first formula is valid outside any 

roles of the two end points are interchanged, and A, containing the positive real axis, and the seco 
is dominant at the left end point If ” is a root of for \are s\<27r/3. In the common part of 
the equation two sectors the two expressions are equly 

In (35) the powers (—s)~* and 2(—s)*“/3 are 
taken positive on the negative real axis and ¢o, 
tinued analy tically into the whole sector of va lit 
In (36) the powers s~*“ and s*” are positive on th 
O(r-*), (29) positive real axis and must be continued analytical! 
into jarg s|<2xr/3. The derivation of (35) and (2¢ 
is omitted except for the remark that actuall, onl 
sign of Im ¢ for large Reynolds numbers to inequali- formula (35) need be proved directly, for fre; 
ties concerning the position of the roots of eq (25) formula (33) one obtains analogous formulas fy, 
and (28). In order to study this question, it is con- gj\) and g3(m), and the identity 


venient to introduce the variable 


0, (28 


then there exist solutions of F=0 of the form 
a 


The preceding results reduce the question of the 
( 


91(n) +92(n) +95(n) =2 x71, 
which follows directly from (23), leads then to 
where . 
; Differentiation of the relation (34) shows that 
and the function 
y’ (s) = Ail 

where Ai(s) is Airy’s integral (ef. [14] and [8] 

5), the function ¥(s) tends to zero as so in t] 


Since the definition (23 implies that (35 
sector 27/3 <arg s<42/3. Integrating (38) one finds 


92(n) =9i (wn) =9s(w"n), (33) therefore, 

all three functions a(n can be simply expressed in y @ | Ai tydt. 
terms of ¥(s). rs 

The values ¢(a@,A) given by (27) have, for suffi- 
ciently large \, positive imaginary parts, if, and 
only if, Re m>0, that is if the zeros of ¥(s) lie in the 
half-plane «/6<args<7/6. Similarly, the values 
¢(a,d) given by (29) have, for large \ positive imagi- 
nary parts, if Re », >0. From (33) it follows that 
the necessary and sufficient condition for this to 
happen is that the zeros of ¢(s) lie in the half-plane 2 

7r/6<arg 8 r/6. Hence all disturbances are ¥(s)=1 3+ | Ai(—t)dt. 
stable , for sufhicre ntly large Re ynolds 7 umbers. if all roots 0 
of J(s) lie in the sector \arg s 7/6. . — a : a. f 
In order to prove that this is indeed the case, a } — Airy’s integral is related to Hankel’s funet 
few properties of the function ¥(s) must first be the identity 


collected A.(s) = \8 Hy, (2! s *), 


From (23), (30), (82) and (33,, it is seen that the ° 3 
“Vy” 


Sed 


Where I, is any path starting at infinity in the seet 
2r/3< argt<4x/3 and ending at t=s. Sine 
Airy integral is real for real s, it follows from 
that ¥(s) 1s also real on the real axis. 

Equations (33) and (37) imply that y(0 
Hence, (39 may be replaced by 


integral representation 
; . where the branch of the three-valued Hankel fu 
Vis)\=—— 93(8)=5 er Sdt (34) tion is determined by the condition that Ai(s) mus 
ont one gol be real on the real axis, the integral in (39) can | 


‘ansformed 1 
is valid. By a calculation very similar to that of “®" formed into 


A;(z,X) in [9] asymptotic expressions for ¥(s) can be w “ aie ) 
derived when |s! is large. The calculation results in vis 9; ./2 aT )EF. 
2143 J 
the following formulas 
Here 


f 


and J; is a path that begins at r=f¢ and goes 
infinity in the upper half plane. As the funetion 
H})(r) has a branch point at sr =0, while p(s) is sing 
valued, the position of Z; with respect to the orig 
must be specified. If 


0O<args<.r, 
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te 


determination of ¢ in (42) is taken for which 


arg [<= (44) 


— 


argument, omitted here, shows that Z; must 
r<arg tT<_2/2. 
» the zeros of ¥(s), it follows from the asymp- 
presentation (35) that outside a sufficiently 
ircle they must all lie arbitrarily close to the 
» real axis. 
he other band, ¥(s) has no real zeros. To show 
ne can use the fact that ¥(s) is a solution of 
ferential equation 


Vv" +sy' =0. (45) 


This can be verified, either by means of (34), or b 


sing (39) and the differential equation for Airy 
«tion. If (45) is multiplied by 2y, the relation 


y 
’s 


1 


‘ l 


_ (OW p- 2 oo 4 _ (y)- . 
ds ait ae ds v ° 46) 


sobtained. Assume, now, that s 
egrating (46) from — 


ais a zero of p(s). 
© toa, it follows then that 


y’"(a)-+4 " ¥ (t)dt=0. 


[his is impossible for real a, because ¥(s) is real on 
he real axis. 
\fter these preparations, it will now be proved that 
he sector #/6<args<11/67 contains no zeros of 
Since ¥(s) assumes conjugate values at conju- 
ate points, only the sector 7/6<arg s<- needs to 
investigated. 
Consider first the portion of this sector for which 


j 


|s| >7o. (47) 
is chosen so large that the remainder denoted 

((s) in formula (35) is numerically less than 
mity, for |s >ro, the absence of zeros of ¥(s) in the 
jomain under consideration is a direct consequence 
{ Rouché’s theorem. An appraisal of the remainder 
based on the following lemma, which is a 

pecial ease of a more general result of Weber [16, 


i, p. 212]. 


LEMMA ] 


Tj R(s) is defined by 


9 


Hi, (8)=(— 


oT 


fo\I 
Ris)|s 


“9 


arg 8s<2r. 


_ Substitution of (48) into (41) followed by an 
integration by parts leads to the expression 


in/4 inet [ a se l : 
: [ _ 1” <2 (5+iR()r ar | 


Or . 
(50) 


v(s)= 5 


From (42), (44) and (47) it follows that ¢ lies in the 
sector —5a/4<arg [<0 of the r-plane at a distance 
from the origin not less than |r|=2r3/7/3. The path 
L; may be chosen as the ray 

t= pe” 
Then 


i, 


s 


th ( aN iR(r)r )dz 


< 


where 


For ro=9, the expression (51) turns out to be less 
than (0.89) |¢"'e% . 


Reinsertion into (50) shows that the quantity O(s 
in (35) is indeed numerically less than one—even 
less than 0.89—for x/6<args<-,\s|>9. There are 
therefore no zeros of ¥(s) in this region. For 
the following extension of Rouché’s theorem is 
needed. 

Lemma 2. Let (2) be regular analytic in the simply 
connected open region G. Let f(x) be defined and 
piecewise continuous on the rectifiable Jordan curve 
Cin G. Assume that on C 


si S70 


¥(2)—f(2)|<\f(z 

Denote by A arg } the total change of arg f(z) as C is 
once described in positive direction. Then the number 
of zeros of W(z) inside C is the integer nearest 
Aarg f/2x. Here arg f(2) can—and must—be defined so 
that it changes continuously at points of continuity of 
f(2) and that it jumps by less than x at the discon- 
tinuities of f(z). 

Proor: The inequality (52) 
bility of ¥(2) having a zero on C. 
by the equation z=A(t), 0<t<1. 
continuous function such that 


excludes the possi- 
Let C be defined 
Define a(t) as a 


a(t)=arg Y(A(t)), O<t< 
Without loss of generality f(A(Q)) can be supposed 


continuous att=0. Because of (52), 


ra 


Re (14+ *—!)>0. 





Hence, 


2) 


7 , 
a())—arg f(A(b) arg ( 1+ 5» (Ost<1) 


determination of 
be called a(t). 
is continuous, 


for a certain uniquely defined 
arg f(h(t) Let this determination 
The last inequality implies that a(¢ 
except for possible jumps of less than x at the dis- 
continuities of f(A) ), and that 


a(t) a(0)] la(t a(0)])< Tr. 


Since the exact number of roots of ¥ inside C is 
a(t) —a(0)]/2x, the lemma is proved 

This lemma will be applied to a contour T in the 
s-plane consisting of three ares T’,, T:, T's as follows: 
’, is the segment from s=0 to |s =9 on the ray 
is the circular are defined by s=9e"*, 
and T; is the segment on the real s-axis 
between 8 9 and s=0. On T,, the function f(s) 
is to be an appropriate partial sum of the ascending 
power series for ¥(s). On T,, the leading term of the 
asymptotic series for ¥(s) can be taken as f(s). On 
';, the function f(s) is identical with y(s). 

The ascending series for ¥(s) follows from (40) 
by termwise integration of the series for Ai(—?), see 
[15] p. B17 Separation of the real and imaginary 
parts for arg s=/6 produces then two series with 
real terms of alternating sign and decreasing moduli, 
at least from a certain term onward. This fact 
makes it easy to find the number of terms needed 
on this ray for f(s) to satisfy the inequality (52). 

The necessary calculations were first carried out 
at the National Physical Laboratory, Teddington 

England) under the direction of F. W. Olver. In 
this calculation ry was taken equal to 10, and 13 
terms of the series were needed to obtain arg f(s) to 


ri: Ty, 


arg 8 


r/i<s0<x, 


two decimals of a radian on T 

Equivalent computations were then performed, by 
a somewhat different method, at the Institute for 
Numerical Analysis of the National Bureau of 
Standards, under the direction of G. Blanch. For 
s'<6 the values of f(s) on T, were obtained by 
numerical integration from tables of the Harvard 
Computation Laboratory [18]. For 6</s\<9 the 
procedure was the same as Olver’s. These computa- 
show that the contribution A, arg f of T, to the total 
change A arg f of the argument of f(s) is —3.42z. 
On his longer ray Olver obtained the value of 

4.11%” radians 

On T, the inequality (52) is satisfied by virtue of 
the appraisal of the remainder used for s >9. The 
contribution A, arg f of T, to A arg f is 


A, arg f 


The corresponding value for |s = 10 is 4.12r. 

Finally, the contribution A; arg f is zero, because 
¥(s) is real on the real axis. The condition (52) is 
here satisfied, since ¥(s) was proved to have no zeros 
on the real axis. 


Summarizing, it is seen that Af=0.01r 
contour formed with the are |s|=9 as well a 
contour with are |s|=10 used by Olver. R 
to lemma 2 then completes the proof that y 
no zeros in the sector x/6 <arg s<11/6r. 

The author expresses his thanks for the ad 4 
help he received from Dr. Olver and from Dr. Bly, 
in this section of the paper. 


4. The Differential Equation 
u™ + ?%z(u"+u)=0 


The occurrence of an inner friction layer js { 
quently regarded as a mathematical consequen 
the fact that the inviscid equation is singular ai 
point where w(y)—ce=0. In the case of plane Coy 
flow the inviscid equation, ¢’’—a’¢=0 is regula 
that point, and it seems plausible that this fy 
responsible for the absence of an inner friction |, 
implied by the unrestricted validity of formulas 
and (11). It will, however, be shown in this sec 
that this view is not quite correct. The total abse; 
of an inner friction layer is actually a very sp 
feature of the differential equation (7). This 
become apparent from the discussion below, of 
differential equation 

u™ + 22(u’" 4 0. 

The reduced equation u’’+-u=0 corresponding 
that differential equation is also regular in the w! 
finite plane. But even a qualitative a priori ar 
ment shows that some form of an inner friction ls 
is to be expected. The solutions of (53) are all reg 
lar at 0, and their fourth derivatives must tly 
fore vanish there. The reduced equation u’ 
has, however, only one linearly independent solu 
with this property. No solution u(z,d) of the | 
equation can, therefore, have a fourth derivative | 
converges to the fourth derivative of, say, cos 
a full neighborhood of z=0. But then w(z,d) car 
even converge uniformly on any circle const 
A>, for u(z,A) is an analytic function of 
the circle and it would follow that wu“ (2, 
uniformly in a neighborhood of z=0, which is a 
tradiction. This remark will now be confirmed 
amplified by actually solving the differential 
tion (53). 

The usual technique of integration by Lap 
contour integrals leads here to solutions of (5 
the form 


u(2,A 


where 


Here € must be a contour such that the intez 
exists, and that the contributions from the bounds 
occurring in the calculations leading to (54) vanis 

The discussion of these paths in this example * 
made complicated by the fact that the integrand 
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+ In order to bring out the 
with the investigations in [9] seven paths 
introduced, which are illustrated in the figures 
t under the assumption that A>0. 
to the determination of the multivalued inte- 
to be taken the following stipulations are made: 
[wo points of the Riemann surface of the inte- 
| are considered to be on the same sheet, if 
can be joined by an are on the Riemann surface 
does not cross the segment connecting the 
-iandt=i. (2) C(V); C(A,), (k 
ie lower edge of C(U;), the upper edge of C(U2), 
nd the right edge of C(U;) lie in the same sheet. 
\s before, the functions defined by (54) will be 
ralled A,, Uy, V according as C is taken to be C(A,), 
,) or C(V), respectively. 
From these definitions the following identities can 


ve derived 


points at t +i. 


oints ft 1,2,3), 


(56) 

(57) 

(58) 

They result directly from an inspection of the paths 
oncerned and the fact that the integrand in (54) 
s nfultiplied by e*”, if the point t=7 is circled once 
hn positive direction. The further identity 
U;— U,=(e"—1) (V—Az) (59) 


ollows by a linear combination of the preceding ones. 
lhe asymptotic calculation of the A,(z,A) for 





For V(z,A) one finds without difficulty the asymp- 
totic expression 


V(2,4 | 1 +#7)— edt 4+-O(r- 
co 


2mri sin 2+O(d~*), 
valid in any bounded domain. 

For the asymptotic calculation of U’, for z in S—S, 
and |t|/>>0, one can proceed essentially as in 


[10], p. 12. Since 


one has 


d™ 


dz” 


— [eres t+arctan f) __ - gf. 61) 
C(U;,) | g2 


The expression in brackets in the last integral is 

numerically less than a constant multiple of 

\~?(/3—t+arctan, t)| max ee h*(3—t+aretany! | hence 
teC(U;) 

the integral is of the form 


. 


‘| t™*1e" f(t, r)dt, 
C(U,) 


where f(t,A) is uniformly bounded on C(L,). Split 
C(U,) into a circle C; of radius r about t=i and two 
rays C, and C, going to, and coming from, infinity 





at the same angle, in two different sheets of the 
The contribution of C; to the 
The contributions of C, and C, 
Hence their sum is 


Riemann surface 
integral is O(A~?). 
differ only by the factor —« 


wih 


>\dt. 


By means of the transformation +=tz, the value of 
this integral is found, similarly as in [10], p. 13, to be 
—0O- O(d? +2 B-4E-2(m42)8) Hence, 


(62) 


for z in S—S, and |f) > >0,m>0. The remainder 
in (62) is O(A~*) for all m, if z is bounded away from 
zero. Otherwise, the remainder tends to zero for 
m <4, but not necessarily for m>4. 

In order to deal with the higher derivatives and 
also with the case |f!<£, observe that the contribu- 
tion to (54) from the path C; is in any case equal to 


a" 


+O (X~*). 
dz™ 


The contribution of the two rays C, and C; is, if ¢ is 
defined by A~?=¢'*, 


yn exp o% 
C3 
+ no} 


o+d~ arctan o ° 


» 


and a short argument, omitted here, shows that this 
Is equal to 


‘da {1 t Or =/3)] 


2(™ 
lo 


where C} is a ray in the o-plane along which Re o°<0 
and whose initial point has at least the distance 
IA?” from the origin. 

In S,, finally, the asymptotic behavior of U,(z,\) 
follows immediately from the identities (57) to (59), 
the asymptotic properties of the A,(2,A), and (62). 
It is then seen that U,(2,A) diverges for || >& in the 
whole interior of S, in the same fashion as the 
corresponding functions in [9] and [10], except for a 


factor e*””’—1=O0(d 


Thus, there is in this problem an 
layer’ 
layer at z=0. 


“nner 
of finite thickness for \t| >&, and also « 


The latter occurs only in the di 


of order four or higher. 
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Diffusion Length of Thermal Neutrons in Water 


James A. DeJuren and Hyman Rosenwasser 


The thermal-neutror 


diffusion length in a hydrogenous medium can he 


simply deter- 


mined by measuring the variation of the thermal-neutron density at large radii from the 


center of a photoneutron source, 


Che necessary requirement is that the relaxation length 


of the source in a medium be less than the diffusion length, so that diffusion dominates at 


large radii. 
was determined by 
was obtained 


1. Introduction 


nnection with the National Bureau of Stand- 
wram of absolute calibration of a standard- 
source, the diffusion length of thermal 

ns in water was measured. A wide spread of 

has been published, and the present experi- 

; s compared with some of the earlier results. 
‘or a high-energy neutron the thermal- 
itron density in water at large radii from the source 
ases in accordance with the relaxation length 

ed below. If a source having only low-energy 
ms is employed so that the rate of neutron 
hermalization per unit volume falls off much more 
with increasing radius than the thermal- 
density from a hypothetical thermal-point 
the thermal-neutron density at large radii 
vary as a function of the diffusion length and 
as in the case of the thermal source). 


source 


apidl 
eutron 


lor 


is only 


- *hotoneutron sources have sufficiently low maximum 


( A . . 
ron energies for this purpose, and two sources 


different neutron-energy maxima were selected 
The general diffusion equation is 


nal aving 
r measurements. 


stal 


dn/dt=q—a-+ DV*n, (1) 
ere n is the thermal-neutron density, q is the rate 

thermal-neutron production per unit volume 
thermal-neutron source-strength density), and @ is 
he absorption rate per unit volume. 

nv) Ye; D=x,,v/3, 
here \, and \,, are the thermal-neutron absorption 
nd transport mean free paths, respectively, for the 
verage Maxwellian velocity, v, and J) is the diffusion 
The steady-state equation is 


n/L?4 with L 


(2) 


q D=0, NAr/3, 
L is the diffusion length. This equation has 
amply treated by Fermi and others in standard 


‘ith water as the medium and at large distances 
the source, one can assume with sufficient 
ey that 

g= Ke /*, (3) 


is a constant, r is the distance from the 
enter, and 6 is the relaxation length. 


Using a Ra-Be photoneutron source, the diffusion length of thermal neutrons 
indium-foil measurements at large radii 
rhis result was confirmed in a second experiment with an Sb"*-Be sources 


763 + 0.015 em 


\ value of 2. 


This relation is very nearly obeyed even for a 
source with a broad spectral energy distribution such 
as is generated by Ra=Be(a,n At large r, 6 is 
usually evaluated by measuring resonance activities 
at energies near thermal with Cd-covered foils. 
Then at large 7 


vi K 
Vn- - ¢ 
i? Dr 
the distribution has spherical symmetry, 
Then 


Since 
n “ur. 


du 
di w ) 


, l 
Ven A 


and 
K 

¢ 
Dr 


du u 


dr? ~ DF 


The solution to the equation satisfying the bound- 
ary condition u->0 as ro is (derivation in appen- 


dix): 
1 [« y - | ; dz ~¢ 7 ; az | 5) 
ay di o/ 2° 


L/L. 


. KL 
ur t 
2D 


where a,=1/6—1/L and am=1/b-4 The expo- 


nential integrals (e~*/r)dr=EI(—ar) may be 
ear 


evaluated from ‘Tables of sine, cosine and expo- 
nential integrals,’ NBS Mathematical Tables MT5 
and 6 (U. S. Government Printing Office, Washing- 
ton, D. C., 1947). 

As r->, the solution may be expressed as 


> P-L 
~ LD? — 5?" 

if L>b, the first term will dominate at sufficiently 
large radii (the necessary requirement for evaluation 
of L). 

Since the second term acts as a_ perturbation, 
measurements with sources having different values of 
b serve as a check on the effect of the perturbation at 
the radii where measurements are possible. 

For a Ra-Be(y.n) source the maximum neutron 
energy emitted is only 0.7 Mev (compared to 13 Mev 
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10 








13 5 


termination of the 


data for the de 


tion length 


] Re 


Fiaurt sonance 


for an a,n source), and the slowing-down length is 
short. Beyond 8 cm from the center of the National 
Bureau of Standards photoneutron standard [1] ' 
the indium resonance activity decreases in accordance 
with eq 3, with a relaxation length 6=2.45 em (fig. 
as 


a 


2. Experimental Arrangement 


The Bureau photoneutron source was placed in the 
center of a 42-in. diameter stainless-steel tank filled 
with distilled water. Four indium foils, 2.0 by 3.4 
cm’ in area and 76 mg/cm’ thick, held in light 
aluminum frames, were suspended by wires from a 
metal framework attached to the top of the tank. 
The foils were centered in the horizontal plane 
passing through the source at positions 90 deg apart. 
A cathetometer mounted horizontally above the tank 
made it possible to position the foils accurately to 
within a fraction of a millimeter at the radii em- 
ployed. The source, a 4-cm-diameter Be sphere 
enclosing a l-curie RaBr, capsule, was pin-mounted 
on a vertical stainless-steel rod extending from the 
tank floor. 

Periodic measurements of the water temperature 
gave a mean value of 22° C, with maximum devia- 
tions of 2 deg. 


3. Counting Techniques 


The foils were activated to saturation at radii 14, 
16, 18, 20, and 22 cm from the origin with and without 
0.020-in. Cd covers. 


! Figures in brackets indicate the literature references 


All four foils were simultane- 


at the end of this paper. 


re ys 


LIXing 


ously counted with “thin” mica end-window 
counters, the outputs of which were fed to a 
circuit, which permitted the foils to be count: 
unit by a single scaler. The sides of the foils {; 
the source also faced the counter windows 
placed in their respective holders inside tl 
Geiger-counter shields. Uranium standard 
used to monitor the Geiger counters, and rey 
cycles were made of standard, background, an 
with counting intervals of 5, 10, and 20 minutes 
respectively, during the decay of the foils. 


wher 


tit 


t 
LOLs 


4. Results 


The data are given in table 1. A correction 
percent must be made for the absorption of 
In resonance neutrons by the Cd covers 
given are determined by the consistency of the d 
as well as the counting statistics. The activi 
given are the counting rates of the foils at saturat 
normalized to the same monitor reading. Fig 
is a semilog plot of radius times thermal activity 


TABLE 1 Data from NBS standard Ra-Be photone r 


Thermal 
activity 
(Au 
Counts 

per minute 


Resonance 
activity 

1.06(Counts 
per minute 


Activity 
without Cd 
(Counts 
per minute 


® Calculated 
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14 





Fiat Indium foils in water. 


Ra-Be source 
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n of r. Beyond 16 em, rA,» decreases 


ntially. 


ToT; 


L , 
« is the experimental error of the ratio 
Air; 
in (“""') te 
¢ the data at 16 and 22 cm as an example, 
5571 
In 636 


€ ) v 
tsi7) ye eo 
ae 
aturat 
. L 


ure 


2.765 + 0.023 cm. 
iVity as 

weighted least-squares fit of the data between 
nd 22 em yield Z=2.763+0.015 cm. 


5. Sb"*-Be Experiment 


\n antimony pellet that was irradiated in the Oak 
dge National Laboratory pile for 2 months was 
substituted for the RaBr, capsule in a 4-cm-diameter 
bervilium sphere. The 60-day Sb™ has three y-rays 
above the Be photoneutron threshold, the principal 
one having an energy of 1.7 Mev and giving rise to 
24-key neutrons [2]. As the neutron energies are 
considerably lower than those from the NBS standard 
moderation occurs closer to the source, and 
relaxation length is shorter. A value of 2.1 cm 
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j 
18 20 





r,cm 


Figure 3. Indium foils in water. 


TARLE 2. 


Data from Sb™-Be photoneutron source 


Thermal 
activity 
fea 
Counts 
per minute 


Resonance 
MIVILY 

1.06(Counts 
per minute 


Activity 
Position |without Cd 
rinem Counts 
per minute 


14.0 
16.0 
18.0 
20.0 


* Calculated, 


for 6 was obtained from the Sb-Be source in water, 
the measurements being made with Cd-covered 
indium and gold foils at radii beyond 8 cm from the 
source (fig. 1,b). 

Indium foils that had been intercalibrated in the 
same slow neutron flux were simultaneously exposed 
to saturation at different radii from the source. The 
foils were successively removed and counted with the 
same end-window tube. For radii greater than 12 
cm, rAy decreases exponentially. At 14 cm the 
corrected cadmium reading is | percent of the thermal 
reading and decreases to 0.35 percent at 20 cm 
(calculated by using b=2.1). The least-squares fit 
of the data (table 2) gives a value of L equal to 
2.76 +0.07 (fig. 3), which is in agreement with the 
Ra-Be result, even though the ratio of indium 
resonance to thermal-neutron counting rates is 
considerably lower for the Sb-Be over the range 
covered. 


6. Discussion and Comparison of Results 


For a purely thermal source located at the origin 


and 


* (r) rédr 
x 6L’. 


| nir)redr 
0 


In analogy, the slowing-down distance, Z,, of 
neutrons from a high-energy source is usually esti- 
mated from the mean-square distance obtained from 
resonance-neutron measurements. r?,, =6L?, and 
the approximate relation 


Pn =P, +6L* 
is used to calculate L’. 

Equation (7) treats the resultant thermal-neutron 
distribution as though the high-energy neutrons 
diffused from the source with a diffusion length equal 
to ¥L?+L*. The shape of the radial thermal and 
resonance-neutron density distributions are actually 
Gaussian near the origin. As a result, the thermal 
and resonance r* values give somewhat overestimated 
diffusion lengths, and the value of LZ obtained by 
this method should represent an upper bound to the 
true figure. 
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Table 3 lists some of the early values obtained with 


RaBe(a n 


sources, using eq 7 


TABLE 3. 


} 


{4 
y+lr 


I 
I 
I 
I 
I 


! 


Since eq 7 involves a small difference between two 
numbers, a 1 percent error in either r* term causes 
about a 3-percent error in L 

Berthelot, Cohen, and Reel [8] have determined LZ 
with a Ra-Be(a,n) source in a water tank by a 
different technique. A thin spherical shell, which 
could be filled alternately with H,O (or H,O+ H,SO,) 
and a CdSO, solution, surrounded the source. The 
difference in the neutron densities outside the sphere 
for the two cases arises from the diffusion of thermal 
neutrons present when the sphere is filled with water. 
Measurements with foils gave a best value of 2.77+ 
0.04 em for the thermal-diffusion length 

Larsson [9] recently has made extensive measure- 
ments of the same type at different CdSO, concen- 
trations. An attempt was made to correct for the 
perturbation in the slow-neutron density in the 
immediate vicinity of the sphere and for the attenu- 
ation of high-energy neutrons by the cadmium. He 
obtained a diffusion length of 3.25+0.13 em. How- 
ever, a different interpretation of the data would give 
a lower value. 

Perhaps the simplest method for measuring L is 
afforded by a nuclear reactor. If a water tank is 
placed on top of a thermal column, the neutron den- 
sity, a few transport mean free paths from the bound- 
ary should decrease as ¢~*'”, providing that the di- 
ameter and height of the tank are much greater than 
the diffusion length. Early reactor experiments 
gave a value of about 2.8 em. Sisk [10] has recently 
obtained a value of 2.67 +0.02 cm at 24° C, using the 
Oak Ridge pile and different boundary conditions. 

Probably most of the variation in the values of LZ 
obtained by the different techniques is due to the 
failure to meet experimentally the rigorous boundary 
conditions required. 

In the present experiment an upper limit to L is 
obtained unless the bracketed term in eq 5 can be de- 
termined. The exponential integrals have been 
evaluated, using L=2.76, and C was roughly deter- 
mined for the radium source by calibrating the foils 
in a standard flux. A perturbation of less than 5 
percent (C>1.3) in L should be caused by neglecting 
the exponential integrals. For the Sb-Be source, 
the exponential integrals are much smaller, and neg- 
lecting them shoufd affect the calculation of L less. 
The Sb-Be(y,n) source was fairly weak, and a stronger 
source would permit measurements beyond 20 em, 
where the exponential integrals would have negligible 
effect. 
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of the diffusion equation. 
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8. Appendix 


Solution of diffusion equation: 


os eS K 


dr i Dr“ 


K 


l L)u "De 


(p 1/L)(p- 


l K 


(p 1/L)(p+1/L) Dr‘ 


Solving by partial fractions: 


l 
(p—WL)( p+ i/Lb 


S,=L/2, 


- l 1 e~r/t 
9 — 
en Le 2p| t1/L p-l 1] r 


r/L+Cyetr/lt+ KL/2D 


, Pen rise-V/k Perl 
enrih dr—er/t 
r 


1/L. 


Let a,=1/b—1/L, a,=1/b+4 


The indefinite integrals may be expressed in terms of defi 
exponential integrals 


e-ar r2 »—ay > 
dr= — dy=— dz: 
r r y ar Zz 


To satisfy boundary conditions at infinity C; 


setting C; = KLC/2D. 


KL a *_-e sah “az 
= -r —_ 1 gtr/L ike ™ 
u=nr ap ° ( f. . dzr+e te - 1 


Since less than! 
are capturec 


To evaluate C, K must be determined. 
percent of the neutrons from either source 
above thermal energies 


Q= arf q(r)rdr -4f — Ke-*/*dr, 
re re 





s the neutron emission rate of the source, 
QO=4abKe->, 

xpression q== Ke-"/* r? does not hold down to the 

the source, so a larger value of ro than the 

just be used, For the Ra-Be(y,n 

be between 4.0 and 


source 
source ro is esti- 
t.5 cm from previous data, 
Qetvb/4ab, an absolute measurement of n would 
1 rough determination of C as all other quantities 
ximately known. 


ential integrals for L=2.76: 


ource 
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